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Abstract 

Subject of our investigations is QCD formulated in terms of physical degrees of freedom. Start- 
ing from the Faddeev-Popov procedure, the canonical formulation of QCD is derived for static 
gauges. Particular emphasis is put on obstructions occurring when implementing gauge con- 
ditions and on the concomitant emergence of compact variables and singular fields. A detailed 
analysis of non-perturbative dynamics associated with such exceptional field configurations within 
Coulomb- and axial gauge is described. We present evidence that compact variables generate 
confinement-like phenomena in both gauges and point out the deficiencies in achieving a satisfac- 
tory non-perturbative treatment concerning all variables. Gauge fixed formulations are shown to 
constitute also a useful framework for phenomenological studies. Phenomenological insights into 
the dynamics of Polyakov loops and monopoles in confined and deconfined phases are presented 
within axial gauge QCD. 



1 Introduction 



The existence of a common formal structure in the theory of the fundamental interactions such 
as QED and QCD reflects the close similarity of the corresponding fundamental processes at high 
momenta. This common formal structure remains manifest in perturbation theory. The high precision 
predictions of pcrturbativc QED have found their counterpart in the quantitative verification of 
perturbative QCD predictions at large momenta. The diversity of low energy phenomena on the 
other hand is a much less obvious consequence of the dynamics of gauge theories. While the low 
energy properties of the weak interaction are understood in terms of the Higgs mechanism as the 
basic non-perturbative element, emergence of the characteristic low energy phenomena of QCD, in 
particular of confinement, remains a central topic in the analysis of gauge theories. Unlike the Higgs 
mechanism, confinement is an intrinsic property of Yang-Mills theories. Identification of the seeds of 
confinement and its consequences is not possible within the realm of perturbation theory. 

Unity of formal structures in the theory of fundamental interactions and diversity of low energy 
phenomena are compatible due to the redundancy of variables - an intrinsic property of gauge the- 
ories. In turn, essential differences in the structure of QED and QCD can be expected to become 
manifest when formulating gauge theories in terms of physical, unconstrained variables. Elimination 
of redundant variables in QCD necessarily introduces non-perturbative elements. In general, the 
process of the gauge fixing cannot be carried out completely In QED it is possible to impose a 
global - for all field configurations valid - gauge condition and thereby to eliminate all the redundant 
variables. It is known [Q that in QCD, a large class of gauge conditions leads to gauge ambigui- 
ties i.e. there exist gauge equivalent field configurations which satisfy the imposed gauge condition. 
Elimination of variables has to respect the complex structure of the manifold of gauge orbits. De- 
pending on the gauge condition, this might be achieved by restricting the domain of definition of 
the physical variables or at the expense of introducing coordinate singularities and thus of including 
singular gauge field configurations in gauge fixed formulations. It is tempting to associate with this 
basic geometrical difference between QED and QCD the characteristic phenomenological differences 
between the strong and electromagnetic forces at low energies. Formation of Gribov horizons 0] or 
condensation of magnetic monopoles Q represent two prominent proposals in which the emergence 
of confinement has been linked directly to obstructions when enforcing a global gauge condition. 

We will discuss in this work progress and problems in the study of gauge fixed formulations 
and present an analysis of non-perturbative properties of QCD in Coulomb- and axial gauge. In 
these gauges it is possible to eliminate completely redundant variables and to formulate thereby the 
theory in terms of unconstrained variables without need for introducing gauge fixing terms or ghosts. 
These "static" gauge choices neither involve the time component of the gauge field nor derivatives 
with respect to time and therefore a standard Hamiltonian description of the dynamics with positive 
norm states can be developed. Necessarily, manifest Lorentz-invariance is lost in these gauges; in 
axial gauge rotational invariance is not manifest either. We start these studies with a description 
of the Faddeev-Popov procedure in the quantization of gauge theories and proceed to the canonical 
formulation with an expression for the Hamiltonian of QCD in static gauges as our main result. Our 
discussion will emphasize from the beginning the role of possible obstructions in carrying out this 
procedure. After implementation of the gauge condition, characteristic properties become manifest 
which distinguish for instance the canonical formulations of QCD and QED. While the corresponding 
Hamiltonians both contain "centrifugal" energies arising, as in mechanical systems, from elimination 
of kinetic terms (electric fields) only in QCD these terms are intrinsically dynamical. With these 
results, a framework will have been established which makes the dynamics accessible to approximative 
treatments. In such formulations no difficult constraints imposed by local gauge invariance such as 
Ward-identities have to be observed. This is essential for our dynamical studies of QCD in Coulomb- 
and axial gauge. In general we concentrate our studies on pure Yang-Mills theories. Inclusion of 
matter fields leaves the structure of the gauge fixed theory essentially unchanged in Coulomb- or 
axial gauge. This is different if the gauge condition depends on matter fields or is exclusively given 
in terms of these, as in the unitary gauge of the Georgi-Glashow model. The unitary gauge condition 
gives rise to obstructions in the gauge fixing procedure closely resembling those encountered in the 
axial gauge - or more generally in diagonalization gauges. We include in this formal section a short 
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discussion of this model. 

Representation of QCD in Coulomb gauge has played a very important role in the development of 
gauge fixed formulations |, ]| 0, §. In QED, the Coulomb- gauge is singled out as the gauge in which 
static charges do not radiate. In QCD this is not the case, gluons and color spin remain coupled 
irrespective of the gauge choice. The focus of our discussion of QCD in Coulomb-gauge will be the 
role of Gribov horizons, i.e. the consequences of restrictions in the range of integration over gauge field 
variables. Gribov horizons are invisible in perturbation theory, they however limit large amplitude 
oscillations. We nonetheless will start this discussion with a perturbative calculation and demonstrate 
how asymptotic freedom will emerge in a formulation where only physical variables are present. It will 
be seen that polarization mechanisms are present which involve the instantaneous Coulomb interaction 
and which are induced by the centrifugal terms in the Hamiltonian. They cannot be associated with 
the coupling to intermediate physical states and therefore do not necessarily generate screening. 
These very same polarization mechanisms will subsequently be shown to be severely affected by the 
presence of a Gribov horizon. We will describe Gribov's approach in dealing approximatively with 
the compactness of the gauge fields and we will discuss the dynamical implications. In particular 
we will study the interaction energy of static color charges and indicate the possible effects of the 
compactness on the structure of the vacuum. 

Formulation of QCD in the axial gauge requires a slightly more complicated setting. For a proper 
definition of this gauge, space has to be compact in at least one direction, otherwise this formulation 
will be plagued by infrared singularities. The axial gauge condition eliminates the component of 
the gauge field corresponding to the compact direction up to the eigenvalues of the Polyakov loops 
winding around the finite spatial extension. In axial gauge we will encounter compact variables and 
singular fields. The restriction in the range of integration affects special variables, the eigenvalues of 
the Polyakov loops. The corresponding Gribov horizon has a simple geometrical origin and is known 
explicitly. The formalism suggests to treat these variables in analogy to quantum mechanical parti- 
cles enclosed in an infinite square well. We will describe this procedure and survey the dynamical 
implications which are far reaching, since these variables, although only a small subset of the dy- 
namical variables, serve as order parameter-fields for the phases of non-Abelian Yang-Mills theories. 
Implementation of the axial gauge condition requires a diagonalization of the Polyakov loops which in 
turn gives rise to singular fields. Their structure will be described and, in absence of a viable scheme 
for treating such variables, arguments will be presented concerning their relevance for the dynamics 
of QCD. 

2 Path Integral and Canonical Formulation 

2.1 Yang-Mills Fields 

Central to the following studies is the discussion of the elimination of the redundant variables. We 
start with a brief description of the Faddeev-Popov procedure and will establish for a certain class 
of gauge conditions the connection to the canonical formalism. In first applications we discuss the 
structure of QCD and QED in the Coulomb-gauge and the canonical structure of the Georgi-Glashow 
model in the unitary gauge. 

We use the following notation and conventions. Gauge fields and their color components are related 

by 

A fl (x)=A«(x)^, 

where the color sum runs over the iV 2 — 1 generators of SXJ(N). The covariant derivative, field strength 
tensor and its color components are defined by 
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Action and Lagrangian are 

S[A] = J d*xC = -\j d A x F^ a F^ = ~\ J (PxtrF^F^. 

As indicated by the notation, in general, color labels will be suppressed. For the canonical formulation 
of QCD one introduces chromo- electric and magnetic fields which, as in QED, are related to the field 
strength tensor by 

E ia (x) = -F 0la (x) , B ia {x) =-~e ijk Fi ka (x) . 
Gauge transformations and transformed fields are written as 

(x) - A™ {x) = n (a,) (a m (x) + 1^ fit {x) t (1) 
where the gauge transformations are parametrized by gauge functions a(x) 



Q (x) — exp {iga (x)} — exp |i<?a a (x) r> 

which for gauge transformations close to the identity yields 

A M (X) -> AW (x) = A M (x) - D,a (x) - | [a, D^a] . (2) 

In path integral quantization of gauge theories redundant variables are eliminated by imposing a 
gauge condition 

f[A] = 0, 

which is supposed to eliminate all the gauge copies of a certain field configuration A. In other words, 
the functional / has to be chosen such that for a given field configuration the equation 

f[AM] =o 

determines uniquely the gauge transformation Q. If successful, the set of all gauge equivalent fields, 
the gauge orbit, is represented by exactly one representative and the resulting generating functional 

Z [J] — f d [A] A f [A] S (/ [A] ) e lS[A]+l I d4xJ " A » (3) 



is given as a sum over these gauge orbits. The Faddeev-Popov determinant 

Sf\A Ml 

A f [A] = | detM | , M= *■ (4) 

a=0 



5a 

M(x,y;a,b)= I d 4 zm> J ' (x, z;a,c) D cb (z)S(z - y) , {x, y; a, b) - " ' ' ^ 1 



ensures that the weight of a gauge orbit in the integral is independent of the representative chosen, 
i.e. independent of the gauge function a(x) (cf. Eq.(g)). If Af [A] vanishes, the gauge condition 
exhibits a quadratic or higher order zero. This implies that at this point in function space, the gauge 
condition is satisfied by at least two gauge equivalent configurations. Vanishing of A/ [A] implies the 
existence of zero modes associated with M 

M X o = (5) 

and therefore the gauge choice is ambiguous. The (connected) spaces of gauge fields which make the 
gauge choice ambiguous 

M H = {A\ detM = 0} 
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are called Gribov horizons. Around Gribov horizons, pairs of infinitesimally close gauge equivalent 
fields exist which satisfy the gauge condition. For a field close to the Gribov horizon we write 

A^ = A^ + 0|t , M[A H ]xo = 0, 

and treat as a small quantity. A Gribov copy 

^'n = + a'^ 

must be connected to A^ by a gauge transformation. The gauge function a must therefore satisfy 
(cf.Eq.§) 

= f[A' = A™ ] - f[A] w -m^a (s) + | [a, 

For this discussion the gauge condition has been assumed to be linear in the gauge fields. With the 
Ansatz for the gauge function 

a = Xo + Sx, Sx < Xo, 
with arbitrarily normalized xo an d $X orthogonal to xo, the above condition can be rewritten as 

m"D IM Sx = -ig([a tl ,xo} + l[xo,D^Xo}) with D = D[A H ]. 

This condition can be solved for Sx provided the right hand side of the equation has a vanishing 
projection on xo- The requirement 

(xo, [<V,Xo]) + {Xo,hxo,Diu,Xo\) = 
determines the normalization of xo- Thus the Gribov copy to has been constructed 

% = a M - D »Xo- 

Perturbation theory in aft yields the eigenvalues A of of the Faddeev-Popov operators corresponding 
to this pair of gauge equivalent fields 

A[A«](xo,Xo) = Uo,m^[AW] Xo ) ~ igm ^ X o, [oW.xo])- 

Using the above normalization condition shows the two eigenvalues to be of opposite sign 

X[A']=-X[A}. 

The spectrum of the Faddeev-Popov operator corresponding to the gauge field beyond the Gribov 
horizon possesses a bound state. 

If on the other hand two gauge fields satisfy the gauge condition and are separated by an infinitesimal 
gauge transformation ea(x), then 

Ma(x) = 

and the two fields are separated by a Gribov horizon. The region beyond the horizon thus contains 
gauge copies of fields inside the horizon. In general one therefore needs additional criteria to select 
exactly one representative of the gauge orbits. The structure of Gribov horizons and of the space of 
fields which contain no Gribov copies depends on the choice of the gauge. At this point we do not 
specify the procedure further but rather associate an infinite potential energy V [A] with every gauge 
copy of a configuration which already has been taken into account, i.e. after gauge fixing, the action 
is supposed to contain implicitly this potential energy 

S[A] -» S[A] -J d 4 xV[A}. (6) 

The above expression for the generating functional can serve as a starting point for deriving the 
Hamiltonian of canonical quantization. To this end we restrict the class of gauge conditions by 



5 



assuming that / is independent of the time component of A and does not contain time derivatives of 
the gauge fields, 

/[A] - 0. (7) 

We now convert the functional integral in Eq.(||) into a phase space functional integral. We write the 
electric field components as 

Ffi = d A? - (DiA ) a 

and introduce auxiliary electric field variables to eliminate terms quadratic in Aq. In the resulting 
expression 

Z[J] = Jd[E] d[Afj] S(f[A]) | det(D^-)| • exp ji J eh (~ (E 2 + B 2 ) - E • (A - DA a ) + J^A, 
Aq can be integrated out 

Z[J]= Jd[E,A]S(f[A]) |det(D-^-)| S(D ■ E - J°) • exp ji ^a; ^-E • A - ~ (E 2 + B 2 ) - JA 

In this procedure, the Gauss law appears as a constraint which will be used to eliminate the electric 
field variables which are conjugate to the gauge fields eliminated by the gauge condition /[A]. To 
identify the appropriate electric field variables we decompose the function space Me of electric fields 
as 

M E =M±®M ll (8) 

with kernel M± and quotient space M\\ 



Mi = < E 



'^E X =0J., M\\=M E /M±. 



SA 

Assuming 6f/5A to be regular (cf. Eq.(||)), i.e. 

detmtm^ 1 ^ 0, 
we can define the projector on M\\ (m = (0, m)) 



P,i = m t - 



1 



m 



'' (mm"f) 
and can represent any E|| 

E||=mV- (9) 

For fixed gauge field A, the functional integral over the electric field variables can - up to a field 
dependent normalization AT [A] - be written in terms of integrations over E^ and the scalar field <p 



d[E] = N[A] J c2[E_lME||] J d[<p]S{Ei\\ -mV). 



With the definition 

' d[A]5{.f[A])N[A] = / d[A x ]N[A], 



the final result for the phase space integral 

Z[J] = fd\E±,A±] -expi /dV[-E ± -A ± -~ (Ei + B[AJ 2 ) 

i f ^ (DE ± - J ))" - J ± A ± } + \nN[A ± }. (10) 
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is obtained. This expression identifies — E^ and Aj_ as conjugate variables and defines the Hamilto- 
nian in the presence of external color currents 

H = y ( i 3 x{i(Ei + B[A ± ] 2 ) + i(mt^i- f (DE ± - J°)) 2 + J ± A ± } + 4 lnA>[A ± ]. (11) 
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This Hamiltonian contains the energy density h(E\ + Bj_) of the physical, unconstrained variables. 
The additional term in the Hamiltonian is the kinetic energy of the eliminated variables. It appears 
in exactly the same way as the centrifugal energy when eliminating angular variables in spherically 
symmetric problems of quantum mechanics. As the centrifugal barrier, this additional term always 
generates repulsion. In the quantum mechanical problem, infinite repulsion is obtained at the center 
where the angles are ill defined and, in gauge theories, where the gauge condition is ambiguous, i.e. 
where with Dm' the Faddeev-Popov operator M possesses one or more zero modes. In both cases, the 
configurations where the centrifugal term diverges describe the limits in the range of definition of the 
corresponding dynamical variables. Integration beyond these limits is effectively eliminated by the 
potential V[A] included in the action (cf. Eq.([|)). In quantum mechanics this procedure corresponds 
to implement the boundary condition at r — by the action of an infinite square well VoO(r). With this 
interpretation, it becomes plausible that phase space path-integral and Hamiltonian do not contain 
the Faddeev-Popov determinant which, in the generating functional, suppresses contributions from 
fields on the Gribov horizon. With the action, also the Hamiltonian contains V[A] and we therefore 
conclude that the Hamiltonian acts in the space of reduced wavefunctionals ^\A\. In particular these 
wave-functionals have to vanish along a Gribov horizon 

m\ AeMH = o, 

which separates infinitesimally close Gribov copies. 

The paradigm in the canonical description of gauge theories is provided by QED in the Coulomb- 
gauge and in a first application of the general procedure we discuss the formal structure of Coulomb- 
gauge QCD and QED. In Coulomb-gauge 

/[A]=divA (12) 
one eliminates the longitudinal gauge fields. With 

y*\3-=6 ab d i 5(x-y), mt=V, A f [A] = | det DV | 

the following expression for the generating functional 

Z[J} = J d [A] | det DV | 6 (div A) e lS[A]+l I d ^ A » . 

is obtained; in this and other gauges discussed below the normalization N[A] (cf. Eq.(|lO|)) is field 
independent and has been dropped. The transition to the canonical formalism is conceptually and 
technically simple for gauge conditions linear in A^ such as (|l2|). In Coulomb-gauge A4± and M\\ 
are the space of transverse and longitudinal vectorfields respectively 

X ± = {E tr } = |e divE = o|, M\\ = {E rotE = o| 

From Eq.(p"l|) the Coulomb-gauge Hamiltonian is easily derived 

H = J rf3 x {iE t 2 r + B[A tr ] 2 -i(DE tr -J )J^AJ^(DE tr -J )+JxAx}. (13) 
The corresponding expression for the Maxwell theory is obtained by replacing D — > V 
H M = J d 3 x u(Vtr + B[A tr ] 2 - + J±A ± j , B = curlA ± . 

In Coulomb-gauge QED, static charges do not couple to the radiation field. This is a very specific 
property of Coulomb-gauge electrodynamics. In axial gauge QED, for instance, static charges couple 
to the radiation field although energy conservation prohibits radiation to result from such a coupling. 
In non-Abelian theories coupling of static color charges to the dynamical degrees of freedom cannot 
be eliminated by the gauge choice. The color spin of static charges always gives rise to a non-trivial 
coupling to the "radiation field". In Coulomb-gauge QED, the l/r is read off directly from the 
Hamiltonian. As in the general case, the Coulomb-energy appears as a centrifugal barrier, however 
it is independent of dynamical degrees of freedom. 
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2.2 Matter Fields 

In this concluding part of our formal considerations we describe the modifications which become 
necessary when matter is coupled to the gauge fields. For static gauge conditions of the type (^) the 
modifications are technical. For the important case of fermions in the fundamental representation as 
described by the additional Lagrangian 

C f = i>(i'fD ll -m)il> 
the modifications can be easily derived with the help of the substitution 

-> - 9i>lvi> 
leading to the following form of the Hamiltonian 

h = y"^{i (E 2 +B[AJ 2 )+ i^ m t_^ (DE± _ j0 _^ t ^ 2 

+i/>W-V - gA)ip + mift (3if) + J j_Aj_}. (14) 
i J 

Since the gauge condition has not been affected by the presence of the fermions, the structure of the 
gauge fixed theory is not significantly altered. Charged matter provides sources of the fields which add 
to the external sources The same remarks apply for the case of the Georgi-Glashow model Q when 
using a static gauge condition. In the Georgi-Glashow model bosons in the adjoint representation of 
SU(2) are coupled to SU(2) Yang-Mills fields. The Lagrangian is 

C GG = ^D fl (j>D» ( f>-V(<t?). 

Under a gauge transformation Q (cf. Eqs.(0J|)) the matter field transforms covariantly 

<f>(x)^<f>W(x) =n(x)(f>(x)& (x) . (15) 

In this case it is of interest to generalize the formalism and allow for dependencies of the gauge 
condition on the matter fields <j> 

/ = /[A,0]. 

We proceed as above from the generating functional 

Z[J}= fd [A, 0] A f [A, <f\ S (/ [A, 0] ) e iS ^ +i I , 



with 

A f [A,q>}= det n ,y J 
da 

We derive the phase space form of Z[J] by introducing fields ir(x) conjugate to the boson fields 4> in 
addition to the electric fields E(x). In this way the time component Aq appears again linearly in the 
exponent of the integrand and after integration yields the Gaufi law constraint 

Z[J]= [d[E,A,K,<j>]A f [A,4,]6(f[A,<f>]) S(T>-E-J Q -ig[ir,<f>}) 

exp ji J d A x - E • A - n<j> - i (E 2 + B 2 + tt 2 + (D0) 2 + V((j) 2 )) - Ja| . 

As in the above case of fundamental fermions, if matter independent gauge conditions are used 

/[A,0 = /[A], 
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only minor changes occur and the result can be read off by inspection (cf. Eq. (|l 
H = y'd 3 a;{i(Ei + B[A ± ] 2 + 7r 2 + (D ( /)) 2 + ^(0 2 )) 

+\ (mt^tDEx-iflM]- J°)\ +J^A ± }. 



For the Georgi-Glashow model the option of matter dependent gauge conditions is relevant for dis- 
playing the physics content of the theory. In particular the transformation property ( |l5| ) suggests to 
use the gauge freedom to diagonalize the matter field. It is convenient to implement this condition 
with the help of an auxiliary field p(x). We write the gauge condition as 

f[A,<f>] =f[c/>} =J>-p T ± + f[A% P (x)>0 (16) 

and integrate in Z[J] over p. The diagonalization condition only would not constitute a complete 
gauge fixing; without /[A 3 ] the gauge condition ( [l6| ) is invariant under Abelian gauge transformations 

This U(l) gauge symmetry has to be removed by the subsidiary gauge condition / [A 3 ]. The Faddeev- 
Popov determinant evaluated for fields which satisfy the gauge condition is easily calculated 

«^L„ - - ,)*) + «l„„^ 

and the following expression for the generating functional is obtained 

Z[J] = J d[A,4>] D[p]A f ~[A 3 )S^f[A 3 ])e iS[A ' pT3/2]+i I dAxJf ' A ^ 

The Faddeev-Popov determinant obviously gives rise to the non-trivial measure 

D[p] = l[p 2 (x)g 2 dp(x). (17) 

X 

The transition to the phase space integral can now be performed with the result 

Z[J] = fD\E,A,ir p ,p]A f [A a ] <5(/[A 3 ]) 6((r>--E) 3 -J°> 3 )-exp{iJd 4 x(-E-A-Tr p p 
-\ [E 2 + B 2 + . 2 + (D p ^) 2 + V{?) + E (( DE )° - JQ -°) 2 ] - JA ) }■ 

^ " a=l,2 

After choosing / the Hamiltonian is easily obtained from this expression. Imposing the (Abelian) 
Coulomb-gauge condition 

/ [A 3 ] = divA 3 

straightforward application of the above procedure yields 

H = f d a x{ l - (Ei + B[A J 2 + ir 2 p + V(p 2 ) + £ ((DE) a - J^f 

J $ ? a=l,2 

+ (D p ^) 2 - (J° + ig\E ± , A ± ]) 3 I (j° + ig\E x , Ai]) 3 ) + J±A ± }. 



The dynamical variables appearing in the Hamiltonian are the pair of (one component) scalar fields 
p, Tip and the pair of color vectorfields Aj_ , Ej_ whose longitudinal color 3-components vanish. These 
expressions illustrate our general remarks concerning the structure of gauge fixed theories. In the 
process of elimination of 2 components of the scalar fields and of the longitudinal color 3-component 
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gauge field two centrifugal terms appear representing the kinetic energy of the eliminated variables. 
While the Abelian gauge fixing via / gives rise to a Coulomb-energy as in QED, implementation of 
the local gauge condition (|l6|) - 5f[<p a ]/5a does not contain derivatives - yields a local form of the 
centrifugal term. This locality makes the Faddeev-Popov determinant factorize into contributions 
from every space-time point and the Gribov horizon consists of those configurations which contain 
matter fields with at least one zero. Here the analogy to quantum mechanics is complete with, for 
given x, the value of the adjoint scalar corresponding to the position of a particle moving in 3-space. 
In agreement with our general discussion, the generating functional contains the non-trivial measure 
( |l7| ) which is absent in Hamiltonian which acts in the space of reduced, "radial" wave-functionals. 



3 QCD in Coulomb- Gauge 

Subject of this section is QCD formulated in terms of transverse gauge fields as physical degrees 
of freedom. In particular we will describe the implications of the compactness of these variables in 
Coulomb-gauge QCD. We follow here closely Gribov's pioneering work in which the limitations 
in the range of the dynamical variables have been suggested as origin of confinement in Yang-Mills 
theories. The following discussion will be carried out within the canonical formulation with the 
expression (^) as our starting point. An object of central interest is the interaction energy of two 
static color charges. In comparison with the corresponding expression in the Maxwell theory, the 
external charges 

J 0a (x)=gp a (x)=gp a 1 (x)+gp%(x) 

are coupled to the the charged gluons by the term DE tr . A corresponding coupling with similar 
consequences appears when a charged matter field is present in QED (cf. Eq.(|l4|)). On the other 
hand the modification of the Coulomb-propagator by the Faddeev-Popov operator is characteristic 
for a non- Abelian gauge theory. 



3.1 Asymptotic Freedom 

We begin our study of the dynamics of the external charges by a perturbative calculation of the 
interaction energy [ flO| |. This discussion not only will prepare the ground for the following non- 
perturbative considerations it will also illustrate how the characteristic antiscreening mechanism 
appears in a formulation in which only physical degrees of freedom and corresponding positive norm 
states are present. For evaluation of the interaction energy in time independent perturbation theory, 
we introduce the total color charge density 

= p a (x) - r bc A t b r (x) E t c r (x). (18) 



and write the Coulomb-energy contribution to the Hamiltonian (13) in a QED like form 

2 E i 2 on g = - 2 P ° S A P * s ' 

The expression for the effective charge density 

9 9 / on x 

A Peff= D~V Pt (20) 



suggests to define the induced color charge density 

Pcff = Pt + And 

with the perturbative expansion 

Pind = -(D-V-A) A+(D 1 v _ A) Pt 



- 9 r bc A b tI ■ V l Pt c + g 2 f abc A b tI ■ V i f cde Ai -V^pt + ... (21) 
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E^ E^ E^ 



Figure 1: Interaction energy of two static charges to one loop. Dashed 
lines indicate the instantaneous Coulomb-propagator, wavy lines transverse 
gluons. 



We now calculate to order g the interaction energy. First order time independent perturbation theory 
yields contributions of arbitrarily high order in the coupling constant. In particular to order g A both 
the first and second order term in the expansion of /9; n( j contribute and the following expression for 
the corresponding contribution to the interaction energy is obtained 



i o r „ i 



A£ (1) = -,<r j /V:,^/<; 

= -^Jp a ^P a - 3. 9 4 r fcc r de < o I y (V i P e ) ■ A d u . I A t 6 r (V i P c ) | o ) conn . 

We have indicated in the formula that only connected interaction terms are retained. These g 2 and g 4 
contributions are given by the first two diagrams of Fig. |l|. The standard expansion of the transverse 
gluon fields reads 

A«(x) = £ Jd\ —Lp e(k, X)A a (k, A)e jk x , (22) 

with the transverse polarization vectors e and the amplitudes A a which in turn are given in terms of 
the creation and annihilation operators 

A a (k, A) = -±= (a a (k, A) + a a (-k, A)t) . (23) 

y/ZUJ k 

With LUk — |k| = k, the gluon propagator in Coulomb gauge is obtained 

(0| AU*) KM) 10)=/ (ga ^ ^ (% - ** - (^) 

and yields after Fourier-transform the following result 

1 »„>(, . My I . t 2 -(q-k)W 



g|k- q| : 



. 9 2 |d 3 fcp?(k)^p«(k)|l 



12iV e g 2 A 2 
48tt 2 11 fc 2 



The logarithmic ultraviolet divergence has been regulated by cutting off the momentum integration 
at A. To order g 4 second order perturbation theory also contributes 

a*» _ J 9 >< o i /„ ^ ^ /„, ^ ft i o ) _ V / A * W £ rfOO ^ ta £ . 

The combined result for Ai?*- 1 ) + AE^ leads to the well known expression for the running coupling 
constant 

3 2 (/i) 



l + i^.9 2 Mlnf 
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Running of the coupling constant is a result of the two competing higher order contributions depicted 
in Fig. yj. The third term yields an increase in the coupling constant with increasing momentum. 
Arising in second order perturbation theory, the sign of this contribution is fixed irrespective of the 
dynamics; more generally, coupling to excited states always provides attraction; from this point of 
view it is clear that with quarks taken into account this attractive contribution gets enhanced. In 
a gauge fixed formulation defined exclusively in terms of physical degrees of freedom this universal 
attraction can be avoided only by processes such as given by the second diagram of Fig. 0. They 
involve the instantaneous interaction and cannot be interpreted as a coupling to physical states. 



3.2 Gauge Fields Close to Gribov Horizons 

In using the standard form of the gluon propagator (|2j) as basis of the perturbative calculation we 
have treated the transverse gauge fields as Gaussian variables. However, as has been pointed out 
by Gribov (|], in Coulomb-gauge and in Lorentz-gauge (as in many other gauges), the restriction 
to one gauge copy leads to a restriction in the domain of the dynamical variables, which invalidates 
application of standard perturbation theory. In various investigations subsequent to Gribov's work, 
the structure of the fundamental domain has been studied. In particular it can be seen by simple 
arguments |[l] that the space of gauge field configurations within the Gribov region - the fundamental 
domain which includes At r = - is bounded in all directions. We start with an arbitrary gauge field 
A tr in the Gribov domain, which necessarily satisfies 



d 6 xto a *(x) {-AS - gf acb jA c tl .(x)V) lu{x)" > 0. 

In absence of zero eigenvalues of the Faddeev-Popov operator — DV inside the Gribov horizon and 
due to the positive definiteness of the Laplace operator the above integral cannot change sign when 
varying the constant 7 in the interval 

< 7 < 1. 

On the other hand for sufficiently large 7 one always can achieve that the field dependent term in 
the Faddeev-Popov operator dominates. With the choice 

iu a (x) = p(x)e ikx v a 

with p(x) > and different from only in a region where none of the gauge field components i 
changes sign. We furthermore assume that |k| is chosen large enough to make the variations of p(x) 
negligible. In this way we always can achieve for k of appropriate orientation and for sufficiently large 
7: 

J d 3 xuj a *(x) (-AS ab - gf acb A c u .(x)V) oj(x) h « 7im( U a V)/ abc k / d 3 xp(x) 2 A c tI (x) < 0. 

Thus for every A tr of the Gribov region the field configuration 7A tr is beyond the Gribov horizon if 
I7I exceeds a certain finite value. Furthermore it can be shown that the Gribov region is convex and 
that the Fourier components A a (k, A) (Eq.(|23|)) of a gauge field in the Gribov region (finite volume 
V) are bounded by IQ 

k, A, a 

These basic insights into the geometry of the Gribov region and further detailed |l3| studies have not 
led to significant progress concerning the dynamics since Gribov's work. We therefore follow closely 
Gribov's original dynamical studies and consider the expectation value of the inverse Faddeev-Popov 
operator which has played such an important role in establishing asymptotic freedom in the Coulomb- 
gauge. We define 

G ^ = <°i (rrV) k i°> = <°i E ^ A(k,A ^ (k ' Atr) io) 
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with eigenvalues A and eigenfunctions x\ of the Faddeev-Popov operator. This quantity is the ghost 
propagator which appears in the functional formulation if the Faddeev-Popov determinant is repre- 
sented by a functional integral over Grassmann variables. After an appropriate ultraviolet regular- 
ization, G(k) becomes infinite only if the Faddeev-Popov operator develops a zero eigenvalue. As 
discussed above this happens whenever the fundamental domain is left and more than one Gribov 
copy is taken into account. In a complete dynamical calculation this is prevented by the action of 
the potential V[A] to be included in the Hamiltonian and which becomes infinite beyond the Gribov 
horizon. We are not able to construct V[A); for a qualitative investigation of the effects of the limited 
domain of the dynamical variables one may require G(k) to be finite and analyze the consequences of 
this requirement. This analysis can be carried out perturbatively by expanding the Faddeev-Popov 
operator around the Laplace operator in the same way as in Eq.(21). Assuming only color singlet 
operators to have non vanishing vacuum expectation values, the two leading terms in the expansion 
of G(k) are (the system is defined in a region of volume V) 



-1 



G(k) aa ' = (0| ( ) \0)5 aa ' - g 2 f abc f cb ' a ' 



\ v^Af;-vj :()) 



^+9 2 (0\ (A t b r .viAt-vJ |0) + . 
Using the expansion (E2J23) of the gluon field, we write by summing ladder graphs 



G(k) = ^(l + a(k,A)) + 



V 



k 2 (l-cr(fc,A)) 



(25) 



with 



a(fc,A)^W^(0|L4 b (q,A)| ; 

b,X J V ' 



10! 



(k - q) 2 



1 - 



(k-q)^ 
k 2 q 2 



We evaluate first a{k = 0, A) in the perturbative vacuum; ground state wave-functional and energy 
are given by 



MA] = ex pj-^7 E ^4 a (k,A)A a (-k,A)j 



with 



k,a,A 



■v 



k, A, a 



d 2 



dA a (k,X)dA a (~k,X) V 



k 2 A a (k, A) A a (-k,X) 



Using the virial theorem, we read off the perturbative result for a(k = 0, A) 



a (0,A) 



11 
3 V 



Y 1 



k,A,a 



Thus in the perturbative vacuum a diverges logarithmically both in the infrared and ultraviolet. 
Independent of the approximations used, it is plausible that only after an ultraviolet regularization, 
a finite result can emerge. Scale invariance of the classical Lagrangian implies that Gribov horizons 
occur on all scales. With Atr(a;) denoting a gauge field on the Gribov horizon, i.e. giving rise to a 
zero mode of the Faddeev-Popov operator xo(x), the rescaled field AAt r (Ax) also possesses a zero 
mode (xo(Ax)) and thus also belongs to the Gribov horizon. For making the above expression finite 
also an infrared cutoff would be required which indicates that in the true vacuum long wavelength 
gauge fields must be suppressed. Within the perturbative approach (Eq.(p5|)) a(k,A) has to satisfy 



a{k,A) < 1. 
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In order to characterize in a semiquantitative way this suppression of long- wavelength excitations we 
modify this condition. As we have seen, perturbative fluctuations favor large values of a and therefore 
one might expect the dominant contributions to arise actually from the region around the endpoint 
and we thus impose the constraint 

o-(0,A) = l. (26) 

As the final result will show, the above constraint for k ^ will be met provided the modified 
constraint (|2^) for k — is satisfied. We treat this condition by introducing a Lagrange multiplier R, 
i.e. the Hamiltonian is given by 

H K = Ho + R^g 2 ]T A a (k, A) A a (-k, A). (27) 



The substitution 



1,2 ]2 , ^_ K 4 _ 4 2 - j_ 



relates H K and Hq] the ground state energy of the constrained system is therefore 



^£V k2+ ^- (28) 

kAa 

The variational parameter R is chosen such that the constraint is satisfied. Using once more the virial 
theorem we calculate 



(0 K |A a (k, A) A a (k, A) | K 
The constraint (Gq) is satisfied provided k is given by 



fc 



2 2 f d 3 k 1 2N c g 2 , 2A 2 

AT ^ I ° J 1^. 



3 J (2-n-)- 3 fcVfc 4 + k 4 37H /€ 

with A regulating the ultraviolet divergence. This schematic calculation explicitly displays suppres- 
sion of the elementary excitations at long wavelengths as a result of the restriction in the range of 
integration of the gauge fields. The energies of the elementary excitations diverge with infinite wave- 
length (Eq.(|28|)). Gaussian fluctuations of massless gluons would extend beyond the Gribov horizon 
and would lead to an overcounting of certain field configurations. This model calculation yields for 



G(k) the result (cf. Eq.(25)) 



k Vvn~i tv A\ AT 2 f A fc 2 ~ 2k ■ q 1 / (k-q) 
-G (k) « 1 - a(k, A) = N c g / — — — — — — - 1 



V w V ' ' ™ J (2tt)3 q ^T^ (k-q) 2 V *V 

with the long wave-length limit 



1 - a(k,A) 



5N c g 2 k 2 lnn 2 /k 2 
k^o 3 k 2 



Thus at infinite wave-length the ghost propagator diverges 

fe^o fc 4 hipr 

while the gluon propagator (Eq.(p4])) vanishes linearly with k 

d 3 TC -*- x (0|yl t o ri (x)4 ri (0)|0) = 1 - ^) <5 afc . 
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It is remarkable that a similar low momentum behavior has recently been obtained from a technically 
sophisticated calculation for Lorentz gauge QCD within a Schwinger-Dyson approach fbi] , |l5| and 
this behavior also seems to be qualitatively confirmed by results from lattice QCD jl6], |l7| . 

The ghost propagator is the essential quantity in the computation of the interaction energy of 
static charges (cf. Eqs.([l9 20)). If we neglect in a first step any polarization of the system when 



introducing a pair of static charges, this interaction energy is, within this model, given by 

V(xi - x 2 ) « J d?k^-^G 2 {k)k 2 



where we have used point charges 



Pl,2(x) = S(x - Xi,2). 



The asymptotics of V(r) at large separations of the static charges is obtained from the long wave- 
length limit of G(k) (Eq.©) 



V{r) 



r 



The limited range in the integration over the gauge fields yields a confining interaction of static 
charges although the increase in energy with distance is too strong. The procedure of introducing 
static charges and neglecting any response of the system is too rough. The polarization of the 
medium should ultimately give rise to the formation of a flux tube and thereby lower the total 
energy of the system. Irrespective of the detailed mechanism, polarization effects will lower the 
exponent in the interaction energy from this value 3. A full calculation is, even within a harmonic 
approximation, technically rather involved. A qualitative idea about the polarization effects might 
be obtained by employing concepts from macroscopic electrostatics. As the perturbative calculation 
indicates, the essential polarization effect - the antiscreening leading to asymptotic freedom - is related 
to the appearance of the Faddeev-Popov operator in Eqs. (|l9| , po| ), only the much smaller screening 
contribution is due to the presence of the gluonic charge in p t (Eq.(|l8|)) which in the following will 
therefore be neglected. Eqs.([l9],|o|) suggest to introduce the Maxwell displacement field Dm with 
Pt as its source while the source of the electric field is p c ff, i.e. E is generated by by external and 
polarization charge. In this interpretation, the limited range of integration gives rise to a dielectric 
constant e which following the standard electrostatic definition 

D M = eE 

is given by 

*>-^- 

Assuming such a linear relation, the energy required to separate two static charges in the polarizable 
medium of the Yang-Mills ground state is 

V(xi - x 2 ) » i Jd 3 xB M -E = J d 3 ke lk ^-^G(k). 

Polarization effects reduce as expected the interaction energy and actually lead to a linear increase 
in the interaction energy 

V(r) « ar 

with the string constant determined by the parameter n (cf. Eq.([29|)) 

3k 2 



807riV c In !f 



The effect of the limited range of integration of the gauge fields due to the presence of a Gribov horizon 
leads naturally to dia-electric behavior of the QCD vacuum over large distances. Such behavior has 
been argued to be responsible for description of confinement in terms of bag- like models Q . Dia- 
electric behavior also implies, by covariance, that the magnetic permeability tends to infinity for 
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small momenta and therefore points to an increasing tendency of the QCD vacuum for spontaneous 
generation of magnetic fields at larger and larger scales. It thus appears that there is a deep connection 
between existence of Gribov horizons and magnetic instability of the QCD vacuum fH| . It is however 
unlikely that, within the above approximations in treating the restrictions in the range of the gauge 
field variables, a qualitatively correct description of the magnetic field dominated QCD vacuum [p0| 
can emerge. We note that in the ground state of the Hamiltonian (p7j), the electric field energy 
dominates, we have 



kAa kAa \I K + ^ 

This model of QCD thus implies a departure from equipartition of electric and magnetic field energy. 
In contrast to the QCD vacuum, this modified vacuum is electric field dominated. This failure is 
most likely due to some inappropriateness in coping with the restrictions in the range of integration 
of the gauge field variables. In our discussion of the axial gauge a different procedure is suggested in 
which compactness is actually the origin of the emergence of a magnetic vacuum. 



4 Axial Gauge QCD 

4.1 Generating Functional and Hamiltonian 

The axial gauge choice singles out a certain direction in space-time. Therefore in this gauge not only 
manifest relativistic covariance is lost, as in any canonical formulation, but in general also manifest 
rotational invariance. Nevertheless, for a variety of investigations and applications, the axial gauge 
representation of QCD is useful. For instance, the physics may single out a specific direction and 
thereby make the axial gauge choice particularly suited such as in the evaluation of interaction 
energy of static charges. The axial gauge choice is the natural choice for finite temperature QCD. In 
light-cone quantization the simplicity of the formalism is intimately related to the light-cone gauge, 
another type of axial gauge. For our discussion it is of particular importance that the gauge fixing 
procedure to the axial gauge including the computation of the Faddeev-Popov determinant can be 
carried out explicitly and in closed form. In this way origin and emergence of compact variables can 
be studied in detail. This in turn makes the physics consequences transparent and leads to various 
applications. In order to properly define the axial gauge one space-time direction has to be chosen 
compact. In this way certain ambiguities and associated infrared singularities of QCD in the axial 
gauge are avoided. Here a spatial direction, the 3-direction, is assumed to be of finite extent L. 
Unlike the temporal gauge (time compact) this choice offers investigations of axial gauge QCD within 
canonical quantization without precluding applications to finite temperature QCD, as will be seen 
later. For this application it is important to assume the compact X3 variables to be elements of a 
circle. We impose periodic boundary conditions for gauge fields and antiperiodic ones for fermion 
fields. In order to be specific we will restrict the discussion to SU(2) QCD; generalization to SU(iV) 
is straightforward. The following gauge fixing procedure is very similar to the one employed in the 
above derivation of the unitary gauge representation of the Georgi-Glashow model. 
We write the axial gauge condition as 

f[A] = A 3 -a 3 ^- + f[A] (31) 

i.e. in the transformation to axial gauge, the 3-component of the gauge field is replaced by a diagonal 
gauge field 03. This field does not depend on X3 and is compact 

a.3 = a 3 (x±), 0<a 3 (x ± )<—. (32) 

gL 

The notation x = (x±,X3) has been used. As in the Georgi-Glashow model a subsidiary condition 
f[A] has to be introduced in order to remove invariance of f[A] under Abelian, X3-independent and 
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in x± local, color-diagonal transformations. This residual redundancy can be eliminated e.g. with the 
following choice 

13 r 3 



f[A}= V / dxsdtAf 

i=l,2 J ° 



(33) 



With the help of the auxiliary field 03 (x±), the generating functional for QCD in axial gauge is 
written as 

Z[J}= I d[a 3 ] I d[A] A f lA]6(A 3 ~a 3 ^)5(f[A]y s W +t f dixJ " A » 



The Faddeev-Popov determinant (cf. Eq. (^)) factorizes 

A/ [A] = \detD 3 \A f [A] 

and det D 3 can be evaluated in closed form. The eigenvalue equation associated with D 3 for gauge 
fields satisfying the axial gauge condition (|3l| ) reads 

[d 3 S ab - ge ab3 a 3 (x ± )] x b n = Kxl 

with eigenfunctions x( x ) which are periodic in x 3 . For a — 3 and for the coupled system of the 
a = 1, 2 components, the eigenvalues are respectively 



A„ — 



2iim 



2iTm 



±iga 3 (x±_), n = 0,±l,±2, 



L L 

and therefore the following value for the ratio of determinants (with zero modes of d 3 removed) 

2 



det D 3 
(det S3) 3 



n=l 



1 - 



/ gLa 3 (x_i_) 
V 2im 



Y[ 71 sin2 gLa 3 (x j_)/2 



(34) 



is obtained. The functional / which fixes a residual U(l) gauge symmetry is assumed to be linear; 
thus the corresponding contribution to the Faddeev-Popov determinant is field independent and we 
can write the generating functional as 



Z[J] 



D\ 



)[a a ] I d [A' ± ] S (/L4]) e is[A' ± ,a 3 ]+iJ *xJ» A ± „ (35) 
where we have separated the measure 

d[A±] = d[A' ± ]d[a 3 ] 

and we have absorbed the Faddeev-Popov determinant into the measure of the auxiliary field 

D [03] = n sin2 {gLa 3 (x x )/2) 9 (o 3 (a:x)) © {2^/gL - a 3 (x ± )) da 3 (x ± ) . (36) 



We conclude this purely formal derivation of the axial gauge representation of QCD with the con- 
struction of the Hamiltonian. The axial gauge condition gives rise to the following decomposition of 
the space of gauge fields (cf. Eq.(p|)) 



M ± = {E tr } = {E I (1 - 6 a<3 )E«(x)T a + d 3 El(x) T3 = 0} 



Mn 



E 



, jf dx 3 E 3 (x) = o| U |e| jf 



dx 3 (diE 2 (x) - d 2 Ei(x)) = 



From Eqs.(|3J,|33|) we obtain 

smA] 



„, . =6 l3 6 ab 5(x-y) + -(l-5 l3 )5 a3 S b3 d l 6(x ± -y ± ) 
SA b t (y) L 
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and using Eqs.(|9|,|ll|) the axial gauge Hamiltonian can be derived 



i(D (E ± + n) j0 )7^p ( D (E-L + tl) - (37) 



P. 

The 2-dimensional longitudinal field appears when eliminating variables by implementing the sub- 
sidiary gauge condition /; it is given by 

1 r L r 3 



rj(x ± ) = ^±y f dx 3 (D x E± - J° ) s 
L Jo 



2 



4.2 Compact Variables 

In this section we discuss the origin of the appearance of compact variables in the elimination of 
redundant variables and indicate some consequences of the presence of such variables. Starting point 
of the above procedure is the definition of gauge condition (|3l|). In order to justify this condition it 
remains to be shown that any arbitrary field configuration can be transformed by a gauge transfor- 
mation into a configuration satisfying this condition. Given an arbitrary field A(x), we pass to the 
axial gauge representation of QCD by applying the gauge fixing transformation 

n(x) = n D (xj_)(P\x ± )) X3/L Pexp^ig dzA 3 (x ± ,z)} (38) 

where 

P(x_l) =Pexp{i.g^ dx 3 A 3 {x)} (39) 

is the (untraced) Polyakov loop and the symbol P denotes path ordering in the integration over 
x 3 . The axial gauge is reached in 3 steps |2]J ■ In the presence of the third factor only, the gauge 
transformation would eliminate A 3 completely. In order to preserve the periodic boundary conditions 
of the gauge fields the second term reintroduces zero mode fields. They in turn are diagonalized by 
Q d which is defined by 

n D (x ± )P{x ± )tf D {x ± ) = e « £ »3(xx)T 3 /2 i ( 40 ) 

This gauge fixing procedure displays a peculiar property of non-Abelian gauge theories, the appear- 
ance of group elements - here the Polyakov loops - in the process of elimination of Lie-algebra elements, 
the 3-component of the gauge field. The only gauge invariant quantities which can be formed out of 
a particular component of the gauge field are the eigenvalues of the Polyakov loop. In QED, it is the 
non-compact field 

QED: — oo < a 3 (x±) < 00 

which is gauge invariant. A division of 03 into a compact and an integer field would be without 
consequences. The integers are the winding numbers which, for fixed x±, are associated with the 
U(l) — ► U(l) mapping defined by the gauge field A 3 (x). These integers cannot be removed by gauge 
transformations. Physically, this reflects the fact that a non-compact field a 3 (x±) must be present, 
otherwise photons propagating in the 1 — 2 plane would not exist in the Maxwell theory. In axial-gauge 
QCD, gluons polarized in the 3-direction can propagate only as composite objects - built up from Aj_ 
and the compact 03. These observations also imply a geometrical interpretation of the Faddeev-Popov 
determinant calculated in Eq.(Q). The resulting non-trivial measure in Eq.(|3^) associated with the 
compact variables is the Haar measure of the group-elements P(x±) and is therefore given (for fixed 
x±) by the the volume element of the first polar angle on S 3 . This interpretation also makes explicit 
the origin of the restriction in the range of integration (cf. Eq.(|32j)) of a 3 . As implied by our general 
discussion, the Haar measure does not appear in the canonical formulation. The Hamiltonian acts 
on reduced wave-functionals which vanish at the boundaries of the domain of definition 

9[A ± ]=0 if a 3 (xL)=0,^. (41) 
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Full and reduced wave-functional are related to each other by 

*[^x]=n j-i ( a 3(^))^x], 

where 

J(a 3 (x±_)) = sin gLa 3 (x±). 

The kinetic energy operator associated with the compact variables acting on the reduced wave- 
functional 

is transformed into 

K a3 V{A ± ] = - V-— - — - J(a 3 (x ± ))— r*[ii], 

t— 1 J{a 3 (x±)) da 3 {x ± ) oa 3 {x ± ) 

X _l 

as is appropriate for the first polar angle on S 3 . Using the full instead of the reduced wave- functional, 
the boundary condition ([ll|) follows from the above form of the kinetic energy operator K a3 . In 
agreement with our general remarks, the last term of the axial gauge Hamiltonian in ( |37| ) is to be 
interpreted as the centrifugal energy which appears in the separation of the angular variables of S 3 . 
Indeed as can be seen from the explicit form of the eigenvalues, D^ 2 develops a second order pole 
when a 3 (x± ) approaches the boundaries or |^ of the domain of definition. 

Here we have used functional methods to derive the gauge fixed formulation of QCD. It is reassuring 
that the results are the same as obtained in canonical quantization Canonically one quantizes 
QCD in the Weyl gauge A$ — and subsequently implements the Gaufi law in order to reach the 
gauge fixed formulation; this procedure actually leads to the Hamiltonian corresponding to the full 
wave functional In turn, using this canonical formulation as starting point, the generating 

functional can be derived and issues concerning the compact variables as discussed above have been 
studied from a different perspective p2| in 1+1 dimensional QCD. 



4.3 Polyakov-Loops and Center Symmetry 

Our discussion has so far concentrated on the origin of compact variables and their impact on the 
structure of the resulting gauge fixed formulation. The compact variables appearing in axial gauge, 
the Polyakov loop variables ([39]) have a significant dynamical role. Their eigenvalues serve as order 
parameters for the confined and deconfined phases of pure Yang-Mills theories |2^, ^J, [2^] . Expec- 
tation value and correlation functions of these variables are related to the self energy of a single 
static quark and the interaction energy of two static quarks respectively and therefore distinguish 
for instance the gluon plasma from the confining phase. The Polyakov loops acquire this dynamical 
significance for symmetry reasons. As will be shown now, the gauge fixing procedure discussed above 
is actually not complete. A discrete residual gauge symmetry, the center symmetry, is still present 
and the Polyakov loop distinguishes the realization of this symmetry. Under gauge transformations 
£l(x), P(x±) transforms as 

P(xx) -> tt(x x ,L)P(x±)& (x_l,0) . 

The coordinates x = (x x , 0) and x = (x±, L) describe identical points and we require the periodicity 
properties imposed on the field strengths not to change under gauge transformations. This is achieved 
if satisfies 

Q(x x ,L) = c n -n(x x ,0) 

with en being an element of the center of the group. Thus gauge transformations can be classified 
according to the value of cq (±1 in SU(2)). Therefore under gauge transformations 

tr(P(x x )) - tr(caP(x±)) ^ ±tr(P(x ± )). 
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A simple example of an SU(2) gauge transformation w_ with c = 



-1 is 



with the arbitrary unit vector tp. Its effect on an arbitrary gauge field is 

A [u>-\ = ^-rxs/L^-i^-rxs/L _ JL^t^. 

9^ 

This representative w_ can be used to generate any other gauge transformation changing the sign of 
tr(P) by multiplication with a strictly periodic (c = 1) but otherwise arbitrary gauge transformation. 
The decomposition of SU(2) gauge transformations into two classes according to c = ±1 implies a 
decomposition of each gauge orbit O into sub-orbits 0± which are characterized by the sign of the 
Polyakov loop at some fixed x° ± 

A(x) e 0± , if ±tr(P(x5.)) > 0. 

Thus strictly speaking, the trace of the Polyakov loop is not a gauge invariant quantity. Only 
|tr(P(xj_)| is invariant under all gauge transformations. Furthermore the spontaneous breakdown of 
the center symmetry in Yang-Mills theory as it supposedly happens in the transition from the confined 
to the quark-gluon plasma phase constitutes a breakdown of the underlying gauge symmetry. This 
implies that the wave functional describing such a state is different for gauge field configurations which 
belong to 0+ and C_ respectively, and which therefore are connected by gauge transformations such 
as cj_ in Eq. (f42|). These symmetry considerations apply equally well when adjoint matter is coupled 
to the gauge fields. They have been applied to 1+1 dimensional QCD with adjoint fermions p(| 
and have shown to be useful in the characterization of the various phases of the Georgi-Glashow 
model [ p6| , ^tJ. On the other hand, with matter in the fundamental representation present, the 
system is not anymore invariant under center symmetry transformations, which change the boundary 
condition of fields carrying fundamental charges. 

In general we have to expect the center symmetry transformations to be present after gauge fixing. 
More precisely, whenever gauge fixing is carried out exactly and with strictly periodic gauge fixing 
transformations (f2, cq = 1) the resulting formalism must contain transformations which change the 
Polyakov loop 

tr(P(z ± )) -> -tr(P(x ± )). 

A discrete residual gauge symmetry is left and each gauge orbit is represented by two gauge field con- 
figurations. The gauge fixing transformation to the axial gauge ([38]) is periodic and therefore residual 
symmetry transformations must be present. It is straightforward to construct these transformations, 
the center reflections. We use the following definition 

Z — i e i " KT ^l 2 & ™ r ^ l L 

These transformations are reflections and change the sign of the Polyakov loop 

Z 2 = l, c z = -l. 

They do not change the gauge condition (|l]) and are therefore symmetry transformations of the 
gauge fixed generating functional (^5|) and Hamiltonian (pTij). 

The effect of Z on arbitrary gauge fields can be simplified by representing the fields in a spherical 
color basis 

%(x) = + iAl(x))e-^ 3 / L (43) 

and by shifting the Polyakov loop variables 

asixx) - -> a 3 {x ± ). (44) 
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As suggested by this definition we will refer to 4> M as charged and to as neutral components of 
the gauge field. Under center reflections, these fields transform as 

Z: a 3 ^-a 3 , A» - -A* , $t , ( M ^ 3). (45) 

and therefore the trace of the Polyakov loop changes sign 

trP(xx) = sin gLa 3 (x±), Z : trP(x±) -»■ -trP(x±). (46) 

The center symmetry transformation Z acts as (Abelian) charge conjugation with the "photons" 
described by the neutral fields A^(x), a 3 (x±). As will be seen shortly these field redefinitions will 



also simplify the description of the dynamics. The phase change in Eq.(43) makes the charged fields 
antiperiodic 

$ ll (x ± ,a?=L) = -$ M (x±,a? = 0). (47) 

If the center symmetry is realized gLa 3 (x±) has to be distributed symmetrically around the origin. 
With the Polyakov loop defined with respect to a spatial compact direction, the center reflections are 
standard symmetry transformations of the canonical theory. In particular, we can associate with Z 
an operator which commutes with the gauge fixed Hamiltonian of Eq. ( |37j ) 

[H, Z] = 0. 

As a consequence, if the center symmetry is realized, the stationary states can be classified as Z-even 
or Z-odd states. 

We conclude this section with some general remarks concerning the relation between covariant 
field theories at finite extension and finite temperature. For the formulation of the center symmetry 
and definition of the Polyakov loops we have considered QCD in a geometry where the system is 
of finite extent L in a spatial direction X3 in order to be able to interpret the center symmetry in 
the canonical formalism. In finite temperature QCD, such a standard interpretation of the center 
reflections as symmetry transformations is not possible and conceptual difficulties may arise [p8f 
concerning for instance the existence of domain walls. However these difficulties can be circumvented. 
By covariance, relativistic field theories at finite (spatial) extension and finite temperature are related 
to each other. By rotational invariance in the Euclidean the value of the partition function of a system 
with finite extension L in 3 direction and (3 in direction is invariant under the exchange of these 
two extensions 

Z(0,L)=Z(L,ff) , 

provided standard boundary conditions in both time and 3 coordinate are imposed on the fields. As 
a consequence energy density and pressure are related by 

e(0,L) = -p(L,0) . (48) 

For a system of non-interacting particles this relation connects energy density or pressure of the 
Stefan-Boltzmann law with the corresponding quantities measured in the Casimir effect. In QCD, 
covariance implies by Eq. (|48| ) that at zero temperature a confinement-deconfinement transition occurs 
when compressing the QCD vacuum (i.e. decreasing L). From lattice gauge calculations p9] it 
can be inferred that this transition occurs at a critical extension L c rs 0.8 fm in the absence of 
quarks and at L c « 1.3 fm when quarks are included. For extensions smaller than L c , the energy 
density and pressure reach values which are typically 80 % of the corresponding "Casimir" energy 
and pressure. When compressing the system beyond the typical length scales of strong interaction 
physics, correlation functions at transverse momenta or energies \p\ <C l/L are dominated by the 
zero "Matsubara wave-numbers" in 3-direction and, as confirmed by lattice QCD calculations p0| , 
are given by the dimensionally reduced QCD2+1. 

4.4 Gauge Ambiguities and Monopoles 

Apart from the discrete center symmetry transformation described by the charge conjugation Z, 
all other symmetries related to the gauge invariance have been used to eliminate A3. In such a 
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case of a global, non-perturbative gauge fixing we have to expect, as argued above, exceptional field 
configurations to emerge. In transforming to the axial gauge, diagonalization of the Polyakov loops 
(Qd in Eq.(|3^)) is the crucial step of the gauge fixing procedure, in which such singular gauge field 
configurations appear. The diagonalization (|4C|) can also be interpreted as a choice of coordinates in 
color space. With Eq. (|4(]), the color 3-direction is chosen to be that of the Polyakov loop at given 
xj_. For studying the emerging singular fields it is convenient to introduce polar (0) and azimuthal 
(ip) angles in color space, 

P(xx) = po(x±) i + ip(x±)r 

with 

p (x±) = cos(^gLa 3 (x±)) 

1 / cos(0(s_l)) e-M*±)M0(?±))\ 
p{xj_)T = sm( — gLa 3 (x±)) 

2 y e itfi ( x ^sm(0 (x_l)) -cos(9(x±)) J 

With this choice of coordinates, the matrix Qd (x±) can be represented as 

/ e Mx ± ) CQS r e t x \ i 2 \ sin (6) ( x \ / 2 ) \ 

V - sin(0 (x±) /2) e -M*x) C os(9(x±) /2) / 

In general, diagonalization or equivalently the choice of coordinates is not everywhere well defined 
and consequently (coordinate-)singularities occur in the associated transformations of the gauge fields. 
Starting from an everywhere-regular gauge field A4 the transformed field 

A'^ (x) = il D (x± ) A M (x) fij, (arx) + (x x ) (49) 

with 

s fl (x±) = n D (x±) —d fi n[ l (x±) (50) 

is in general singular with 11d- While the homogeneous term can at most be discontinuous the 
inhomogeneous term diverges. We note that s M (x±_) is determined exclusively by the Polyakov- loop 
variables po , p and furthermore it can be shown that the parameters characterizing the singularities 
are given in terms of the gauge invariant eigenvalues of the Polyakov loops po. We represent the 
inhomogeneous term s M in a spherical color basis (cf. Eq.(|43|)) 

s M (x±) = a* (x±) t 3 + -j= (0* (x±) (n - ir 2 ) + h.c), 



where 



Wx±) = =-e- i *(Bm0d M <p-iB ll 6) = — — e~* {$> - iff) ■ d» p (xx) 

2g 2gsmgLa 3 

« r s 1 /-, n \ r, I+COS0 1 _ . 

a^xx) = --(1 + 0080)^ = - y^p(^), (51) 



with the standard choice of unit vectors (in color space) 




cos cos ip 

Cp = I cos ip I and = | cos sin ip 

— sin 9 

These expressions display the nature of the singular fields and describe the conditions under which 
singularities occur in the process of gauge fixing. We assume the gauge field (x) and in particular 
the corresponding Polyakov loop P = (poiP) to be smooth before diagonalization. Singularities 
(poles) occur at points x^' S , where the Polyakov loop passes through the center of the group and 



does not define a direction in color space, 



P(a£' s ) = ±11. (52) 
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This happens if 

gLa 3 (x±) =0 or gLa 3 (x S ± ) = 2n 



i.e. if the Polyakov loop variable reaches the border of the fundamental domain (Eq.(32)). This 
requirement determines a point on the group manifold S 3 and thus, for generic cases, fixes (locally) 

N S 

uniquely the position x ± ' . In 4-space the transformed gauge fields are thus singular on straight lines 
parallel to the 3-axis. As a consequence of the particular gauge condition involving the diagonalization 
of the X3 independent Polyakov loop, they are also independent of X3 and have a vanishing /i = 3 
component. The singularities can be classified according to the value of the Polyakov loop and we 
shall refer to them as north and south pole singularities according to the respective positions of the 
Polyakov loop on the group manifold SU(2) (= S 3 ). Thus we can assign a "north-south" quantum 
number or charge z to this singularity (i.e. the range of z is the center of the group). In addition to 
poles, the field also exhibits (static) string like singularities along the line 9(x±) = representing 
a surface in 4-space. The charged gluon fields too, have poles at x , ' and discontinuities along the 
strings 9(x±) = 0. Although the points x^_' S are characterized by Eq.([52|) in a gauge invariant way 
(by the degeneracy of two gauge invariant eigenvalues), in general those points have no particular 
significance in other gauges. To illustrate this general discussion we consider a particularly simple 
singular field which arises if we identify color orientation and spatial orientations. For this we assume 
the x± space to be Euclidean (imaginary time) and we use a vector notation (e.g. a = (00,01,02)). 
Identification of color and spatial orientation 

y t 

<P (x±) = (fio = arctan — , 9 (x±) = 9q = arccos — . 

x y/ x 2 +y 2 + t 2 

generates a singular field a^(x±) in Eq.(pl|) whose neutral component is the vector potential of a 
Dirac monopole |3lJ of charge 2ir/g, 

i \ 1 l + cos6» . , . 

a(iB_L) = -- ^-7, Vo 53 ) 

2g r sin 9o 

with Abelian (neutral) magnetic field 

b = curl a = — . 

2g x\ 

Associated with the singular neutral component is a singular charged component which is given by 

cj ) ( x± ) = -^-(^ +ie )e- i ^. (54) 

This example exhibits a general property of the singular field configurations generated by the 
diagonalization. The strength of the singularity of the neutral component is determined by the 
winding of the color orientation when the orientation in x± space is appropriately varied. The 
magnetic charge m is quantized and given by 

2-rtn 

m= . n = ±l,±2, ... . 

9 

On the other hand, the structure of the singular charged component is, in general, not determined 
exclusively by topological properties. However the singularities of neutral and charged components 
are intimately related. The Abelian field strength corresponding to the neutral singular fields, a 
central quantity of the so called Abelian projection 

f»u = df.al - <9„a^ , fi, v ^ 3 

is singular at the position of the monopoles. On the other hand, the complete non-Abelian field 
strength built from the inhomogeneous term s\ of Eq.d5G) actually vanishes, 



[s] = d^Sy - d v Sfj, + ig [s M , s„] = 0, 
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i.e. the singular Abelian field strength is exactly canceled by the non-Abelian contribution to F 
generated by the singularities in the charged gluon fields. Thus "Abelian" monopoles have finite or 
possibly even vanishing field strength. Singularities in the gauge fields necessarily cancel in gauge 
invariant quantities; they have been produced as coordinate singularities since we insisted on fixing the 
gauge globally. Such cancellations can be achieved only if the connection between neutral and charged 
singular components of the gauge fields contained in the above expressions ( |49| , |5l] ) is not disturbed. 
Cancellation of the singularities in the Abelian field strength by the non-Abelian commutator must 
happen quite generally, the gauge fixing procedure cannot affect gauge invariant quantities even at 
the positions of the monopoles and along the strings. The above considerations also make clear that, 
in general, there are no bounds on the action associated with singular fields. The simple example of 
the gauge field 

Ap(x) = — t 3 6^3 

illustrates this point. The corresponding action is zero and the whole space will be filled with 
monopoles after transforming to axial gauge. On the other hand, fluctuations around singular field 
configurations will disturb the delicate balance between Abelian and non-Abelian contribution to 
the field strength and, in general, an infinite action will result. Fluctuations have to satisfy specific 
requirements to yield finite action. The charged components of the fluctuations 5(f) have to satisfy 
the condition 

8(j>(x) e 2lv<yX± -^ continuous along the strings (55) 



and both neutral and charged fluctuations have to vanish at the position of the pole. 

Unlike the points Xj_' which are determined in a gauge invariant way (cf. Eq.(p2|)) the location of 
the strings is to a large extent arbitrary and depends on the details of the subsidiary gauge condition 



f[A] of Eq.(3f ). An Abelian rotation with an rr^-dependent gauge function only affects the f[A] part 
of the gauge condition. This subsidiary gauge condition can be used to simplify the description of 
the strings. 

For the following discussion and for further applications it is convenient to make explicit in the 
generating functional the contributions from singular field configurations. This is not necessary if 
a complete evaluation of the path integral is attempted. For approximative evaluation however a 
decomposition into singular fields and fluctuations is useful. To this end, we classify the gauge fields 
according to the number n = (riff ,ris) of north (riff) and south (ns) pole singularities i.e. according 
to the number of times the Polyakov loop passes through the "north" and "south" center respectively 
and denote the corresponding fields by a" . On the basis of this classification the generating functional 
can be decomposed as 

Z[J] =J2z n [J\ =J2f D [^W±] S (f[ A ]) e lS l A '^ a "-* /9L } +l J d4xJ " A± - (56) 

n n J 

As in Eq.(|35|) the integration variables, the unconstrained degrees of freedom, are the 3 components 
of the gauge field A'j_ = ^ 3} and the eigenvalues of the Polyakov loops. For n ^ 

singularities show up in the gauge field components A' ± which in turn are determined by the x± 
dependence of a". One therefore can split A' ± into singular and fluctuating fields 

A' ll (x)=a»(x ± ) + 8A' IA (x). (57) 

For given a 3 (x±), the singular fields 



can be constructed as generalization of the Dirac monopole solutions (Eqs.( |53| , |54| )). For completeness 
we also write down the measure (cf. Eq. (p6|) ) of the Polyakov loop variables, which after the field 
redefinition (El) reads 



D [a 3 ] = J] cos 2 (gLa 3 (x ± )/2) O - <%(x±) ) da 3 (x±_) . (58) 
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4.5 Phenomenology of Polyakov-Loops 

Up to this point our discussion has focused on the development of the formalism for gauge fixed 
theories and a description of the properties of the physical variables reached after gauge fixing. Within 
axial gauge QCD, it has been possible to carry out a rather complete analysis of the properties of 
the physical, unconstraint variables. We have demonstrated the emergence of the compact Polyakov 
loop variables and of the singular field configurations which appear whenever the compact variables 
reach the border of the fundamental domain. Compactness of some of the variables and the presence 
of singular gauge fields not suppressed by an infinite action, constitute characteristic properties of 
the gauge fixed non-Abelian theory. 

In the following sections we will discuss to what extent the non-trivial phenomena of QCD may 
be associated with these properties of the gauge fixed theory. In this endeavor the Polyakov loop 
variables will play a central role. It will be an important asset that, in axial gauge, these variables 
which serve as order parameters occur as elementary rather than composite fields. Therefore in the 
absence of a viable approximation scheme, the known properties of QCD can in turn be used to 
deduce dynamical properties of these variables. 

The central phenomenon of Yang-Mills theories is confinement. Polyakov loops are objects whose 
dynamics is intimately linked to this phenomenon. The spectrum of the Polyakov loops - with a 
spatial direction compact - reflects directly presence or absence of confinement. We consider the 
correlation function of two Polyakov loops separated in Euclidean time. For large separation, this 
correlation function is given by 

lim (0|trP(i E )trP(0)|0) = ce~ E - tE , 

t — >oc 

where, after Wick rotation, we have made the following choice of coordinates 

x± = (i B ,0,0). 

P_ denotes the energy of the lowest state, which can be reached in applying the Polyakov loop 
operator with negative Z-parity to the ground state. On the other hand, after a further rotation 
in the Euclidean, time and 3-axis can, up to a sign, be interchanged. In this operation the value 
of the correlation function does not change; it however now acquires a different interpretation. The 
correlator 

lim (0|trP(a; = t B )trP(0)|0) = e - LV( - x) 

t — >oo 

describes the same system at finite temperature with T = 1/L and the Polyakov loop correlator 
(corresponding to compact Euclidean time), as is well known p4j, is given by the free energy of a 
pair of static charges. In the confined phase 

1/L = T < T c , V(r) = or 

the interaction energy is linearly rising with the slope given by the string constant and therefore we 
conclude 

P_ = oL, if L > L c = 1/T C . (59) 

The lowest energy of states which can be excited by the Polyakov loop operator therefore increases 
linearly with the extension of the system. In particular, the ground state cannot contribute which is 
guaranteed provided the ground state is symmetric under center reflections (cf. Eq.(|46|)) 

(0|P(s_l)|0> = 0, ifZ|0) = |0). 

In the deconfined phase we expect Debye screening to take place giving rise to an interaction energy 
of static charges 

V D ~ \e- mor . 

Following the same line of arguments this implies that, for extensions below the critical value, the 
ground state breaks the center symmetry and the excited states which can be reached by P(x_i_) 
exhibit a gap P_ 

L<L C : Z\0) ^ |0) P_ = m D . 
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A further characterization of the confined phase and the dynamics of the variables 0,3 can be obtained 
through a discussion of adjoint Polyakov loops. The adjoint Polyakov loop is defined with the matrices 
T a of the adjoint representation as 



P ad (z ± ) = ^trPexp \ig jf dzA&x±, z)T a J . 

After gauge fixing also the adjoint loops are given in terms of the variables C13 (cf. Eqs.( |3l]j4^ )) 

P ad (x±) = ^{1 -2 cos gLa 3 {x ± j) . (60) 

Adjoint charges can be screened by gluons; thus for sufficiently large distances, the interaction energy 
must tend exponentially to a constant with the exponential slope determined by the lowest glueball 
mass m g b- This implies a non- vanishing ground state expectation value 

(0|(f-2cos 5j La 3 (x x ))|0) ^0 

and a gap E+ of states which which can be excited by P a <i 

E+ = m gb . (61) 

Being not forbidden by symmetry requirements, a non- vanishing vacuum expectation value is natural. 
These facts concerning the dynamics of the variables as(x±) strongly suggest the following properties 
of the spectrum of states which can be reached by applying operators built from these fields. In 
the confined phase, the ground state of the system is even under center reflections. The spectra of 
excited states depend crucially on the Z-parity. States with positive Z-parity are physical, "hadronic" 
states, i.e. in pure Yang Mills theories they describe e.g. glueballs (with vanishing 3 component of 
the momentum ^3). States with negative Z-parity and vanishing p% exhibit a gap which becomes 
infinite with infinite extension of the system. In other words in this limit these states are frozen and 
we may expect this to be the case also for states with p% ^ but significantly smaller than l/L c . In 
the confined phase no gauge-invariant operator can connect states belonging to the 2 different sectors 
due to their different Z-parity. The presence of quarks will substantially change this picture. With 
decreasing extension, the excitation energy of the Z-odd states decreases and when approaching the 
critical extension at which the confinement-deconfinement phase transition takes place, the gap in 
this sector is of the order of the glueball mass, which at infinite extension characterizes the gap in 
the Z-even sector. More quantitatively we know from lattice calculations [[32| |3^] that in SU(3) with 
T c w 220 MeV, the values of the lowest glueball mass and the gap in the Z-odd sector at T w 170 MeV 
are 

m gb (170 MeV) « 770 MeV, cr(170 MeV)/T c « 670 MeV. 

Thus at the phase transition the continuously decreasing gap suddenly vanishes together with the 
string tension and the Z-odd states become available and contribute to the thermodynamic quantities 
such as pressure or energy density. At the same time the glueballs in the Z-even sectors disappear. 
With the ground state breaking the center symmetry, the two classes of states are now coupled and 
are therefore in thermodynamic equilibrium. Thus the confinement-deconfinement phase transition 
is not just a melting of the glueballs. Rather at the transition a whole sector of the Hilbert space, 
completely decoupled below the phase transition and not accessible to any physical observable, joins 
the physical states in the center-symmetry breaking plasma. 



4.6 Theoretical Approaches to Polyakov Loop Dynamics 

After having described the most prominent properties of the dynamics of the Polyakov loop variables 
we now turn to attempts to provide theoretical understanding of some of these gross features of 
QCD. On the basis of the expression ( pff ) for the generating functional we will describe a hierarchy 
of approximations in the evaluation of Z with increasing complexity. It is clear from the outset that, 
even for modest success, certain non-perturbative elements have to be incorporated. 
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The QCD generating functional in the naive axial (or temporal) gauge is obtained if only the 
sector without singularities is kept and the dependence on the eigenvalues of the Polyakov loops is 
disregarded. As a consequence of these approximations, the generating functional becomes actually 
ill-defined as has been noticed early by Schwinger f34|| . Definition of propagators requires certain "ie" 
prescriptions. In the course of the approximations, the center symmetry got lost. 

Still, keeping the zero singularity sector only one might proceed by accounting for the dependence 
of Z on a 3 . The simplest form of these dynamics results, if these variables are treated as Gaussian 
variables, i.e. if the non-flat measure of Eq.(|5q) is replaced by the flat one 



D [o 3 ] s» \\_daz (a; j.) . 



In this way, one effectively treats the Polyakov loop eigenvalues as the zero modes in QED. It is 
therefore not surprising that the center-symmetry is lost again and and a perturbative picture emerges 
with the phenomenon of Debye screening as the leading dynamical correction to the description of 
QCD as a system of non-interacting gluons ]3q |. 

As we now will discuss in more detail, first characteristic properties of QCD are encountered if, 
still disregarding singular field configurations, the non-flat measure of the Polyakov loop variables is 
properly taken into account. In particular, the perturbative phase reached in this way will be seen 
to be center-symmetric. In the last section we will address the possible role of the singular field 
configurations in axial gauge QCD. 

A crucial element in the following discussion will be the compact, i.e. non-Gaussian nature of the 
Polyakov loop variables a 3 (x±). We have seen that the appearance of compact variables is common to 
most of the formulations of QCD in terms of unconstrained variables. To study the consequences we 
use the canonical formulation and disregard in a first step the coupling of Polyakov loop variables to 
the other degrees of freedom |36|. In the absence of such couplings, the Hamiltonian of the Polyakov 



loop variables reads (cf. Eq.(37)) 



d 2 x 



' 2 ' £ (V fl3 (x.)) 2 



2L<5a 3 (xjJ 2 2 
and if space is discretized (lattice spacing i) 

H =-WT. + -h E (S.C*, + Q - ^(x,)) 2 = V + h" , (62) 

where S — (£, 0), (0, £) denote the fundamental vectors of the lattice and the dynamical variables have 
been rescaled 

a 3 (x±) = gLa 3 (x±)/2 . (63) 

For weak coupling to other degrees of freedom, the Hamiltonian ( j6"2| ) describes both, photons in QED 
and Polyakov loop variables in QCD. It however acts on wave functions belonging to different spaces. 
In QCD the compact nature imposes the boundary condition (cf. Eq.( ]4l| , f44| )) 

A 7T 

^[5,3} = whenever 03 (xj_) = ±— for some . (64) 
In order to display the non-trivial dynamics described by the Hamiltonian h in conjunction with the 



constraint of Eq.(64) we consider first the case of electrodynamics where no such constraint is present. 
As is well known, by discrete Fourier transformation, the elementary excitation can be determined 
and the following dispersion relation for the (lattice) photons, 



a 4 ^ . 2 (SkY , 2tt 



w k = ^2^sin^j— k=— (ni,ra 2 ) (65) 
s 



is obtained, with the standard continuum limit 



2 1 2 
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The ground state wavefunctional 



*[a 3 ] = exp<^ -- / d 2 2 J _d 2 2/_ L a 3 (x_ L )i ; s:(x ± ,x'_ L )a 3 (x'_ L 



(66) 



expresses by its non-locality 



strong correlations in the system. While electric and magnetic field energy contribute equally in the 
normal modes of QED, in QCD as a consequence of the boundary condition ( |3^ ) 

(a 3 (x_L +5)- a 3 (x±)) 2 < tt 2 

the magnetic field energy h m becomes negligible in the continuum limit (L/£ — > oo). Comparable 
contributions from electric and magnetic field energy could result in the continuum limit only if we 
assume linear instead of logarithmic running of the coupling constant 

g 2 ~ £, £^Q. 

Dominance of the electric field energy yields the reduced ground state wavefunctional 

1/2 



cos(a 3 (x_L)) 



(67) 



States of lowest excitation energy are obtained by exciting a degree of freedom at one site xj_ into its 
first excited state; this is achieved by replacing in the ground state wave functional 



and an excitation energy 



cos (a 3 (xjj) -> sin(2a 3 (xx)) 



AE = 



8 £ 2 



(68) 



results. Such excited states can occur at any site and therefore these states are highly degenerate. 
This degeneracy is lifted by the magnetic coupling. A perturbative evaluation yields a "band" of of 
excited states characterized by the discrete momenta (cf. (pq)) and with excitation energies 



A£ k = 



3 g 2 L 64 1 



8 P 9tt 2 g 2 L 



- 1 



The qualitative differences in the structure of the ground states of QED and QCD respectively 
(Eqs.( p6p7| )) imply very different properties of the corresponding elementary excitations. Built on 
the highly correlated QED ground state the photons appear as collective excitations with excitation 
energies vanishing in the long-wavelength limit. In QCD, the elementary excitations are localized in 
configuration space and are due to formation of non- vanishing electric flux. In the absence of couplings 
to the other degrees of freedom, the chromoelectric fields formed in the elementary excitations with 
lowest excitation energy are located on just one transverse lattice site and the flux tube is infinitely 
thin. It winds around the compact 3 direction and thereby gives rise to an excitation energy which 
increases linearly with the extension i.e. in this limit 

E- = (JstcL, 

with the string tension given by the strong coupling lattice result p7f 



M 2 ' 



The ground state of the system is, in the strong coupling limit, an eigenstate of the electric field op- 
erator. In QED with photons propagating in the 1 — 2 plane described by Gaussian variables a 3 (xj_) 
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such states are not normalizable and would entail infinitely large fluctuations in the magnetic field 
energy. Thus the structure of the vacuum concerning these ^-independent fields is very different in 
the Abelian and non-Abelian theory. In QED the virial theorem yields the standard result: mag- 
netic and electric fields contribute equally to each normal mode. In axial gauge QCD the Jacobian 
invalidates equipartition. Chromoelectric a^-independent fields are absent; the square of the reduced 
ground state wave functional ( |67| ) is nothing else than the Jacobian (^8|) and thus the corresponding 
full wave- functional is constant. In turn, the fluctuations in the magnetic field at different lattice sites 
are not correlated and the ground state energy is due exclusively to these uncorrelated magnetic field 
fluctuations. Therefore the corresponding contributions to the "gluon condensate" are dominated by 
the magnetic field 

(5(0|E 2 -B 2 |0) = -£(0|B 2 |0) . 

With the electric flux quantized and resulting in the gap ( |(38| ) in the spectrum, this magnetic field 
dominated ground state of the Polyakov loop variables exhibits a dual Meissner effect, i.e. it resists 
penetration of x 3 independent electric fields pointing in the spatial 3 direction into the medium. 
This is not a result of a condensation of monopoles but a consequence of the compactness of the 
relevant gauge fields. It is interesting that even in the crude approximation of keeping only one 
particular kind of gluonic degrees of freedom, this model displays features which are reminiscent 
of the phenomenology of the "magnetic QCD vacuum" ]2C| ]. The compact nature of the Polyakov 
loop variables appears to reproduce the most striking feature of the spectrum in the Z-parity odd 
sector of the theory (cf. Eq.(p8|)), the linear increase of the gap with the extension. It thus displays 
certain confinement like properties, it gives rise to a magnetic vacuum and exhibits the dual Meissner 
effect. Quantitatively, the description of these variables as infinitely thin flux tubes and the deduced 
strong coupling value of the string constant is unrealistic; accounting for the coupling to the other 
degrees of freedom might be expected to improve these results. On the other hand at this level of 
the development, the formalism fails in not displaying any significant differences in the dynamics of 
confined Z-odd and "hadronic" Z-even states. 



4.7 Perturbative Results for the Confined Phase 

The following studies concerning the effects of the coupling of the Polyakov loops to the other gluonic 
variables will be performed on the basis of the generating functional (|56|). In particular we will 
be interested in the properties of the Polyakov loop correlation functions. For their evaluation, we 
observe that the (in the continuum limit) infinitely large gap (cf. Eq.(|o8|)) prohibits the Polyakov loop 
variables to propagate. Indeed our above results are easily translated into properties of the relevant 
correlation functions. In the absence of coupling to the other degrees of freedom the generating 
functional written for discretized space-time and in terms of the rescaled variables (|63|), is in the 
Euclidean, given by 



Z = J d[a 3 ]exp|-l/2 J d 4 x(d,Mx±)) 2 

= / Ylda 3 (x±) cos 2 a 3 (xjj expj^— (a 3 (y± + S±) - d 3 (y ± )) 2 } . 



In the continuum limit 



v±,s± 



1 £ i L n 

?Z^I ln 7^ ' 



and therefore the nearest neighbor interaction generated by the Abelian field energy of the Polyakov 
loop variables is negligible. As a consequence, in the absence of coupling to other degrees of freedom, 
Polyakov loops are ultralocal, they do not propagate, 



(0|r(o3 (xx)a 3 (0))|fi) ~ (^Y ±/e ^ t 3 (x±) 
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Propagation of excitations induced by a 3 (x±) can only arise by coupling to the other microscopic 
degrees of freedom. Ultralocality permits the Polyakov loop variables 03 to be integrated out and the 
following effective action results 

S cS [Afj] = S YM [A„, ^1 3 = 0] + Sgf [J dz 4] + M 2 j d A x <S>l{x)^{x). (69) 

The Polyakov loop variables have left their signature in the geometrical mass term of the charged 
gluons (cf. Eq.@) 

M 2 — (7r 2 /3 — 2)/L 2 (70) 



and in the antiperiodic boundary conditions (Eq. (|47|)). The neutral gluons remain massless and 
periodic. The antiperiodic boundary conditions reflect the mean value of the Polyakov loop variables, 
the geometrical mass their fluctuations; notice that in both of these corrections, the coupling constant 
has dropped out. We emphasize that periodic boundary conditions for the gluon fields are imposed 
in the original expression (^5|) of the generating functional. The antiperiodic boundary conditions 
implicit to the definition of the generating functional in (|5^) account for the appearance of Aharonov- 
Bohm fluxes in the elimination of the Polyakov loop variables. Periodic charged gluon fields may be 
continued to be used if the differential operator $3 is replaced by 

ds^d 3 + ^[rs, . (71) 

As for a quantum mechanical particle on a circle, such a magnetic flux is technically most easily 
accounted for by an appropriate change in boundary conditions - without changing the original 
periodicity requirements. With regard to the rather unexpected physical consequences, the space- 
time independence of this flux is important, since it induces global changes in the theory. These 
global changes are missed if Polyakov loops are treated as Gaussian variables. 

The role of the order parameter is taken over by the neutral color current in 3-direction u (x± ) 
which is generated by the 3-gluon interaction 

u (x x ) = i [ dx 3 $t (z) g 3 $m (x ) . ( 72 ) 
Jo 

This composite field is odd under center reflections (cf. (E5h) 

Z : u(x±) —>■ —u(xj_). 
It determines the vacuum expectation value of the Polyakov loops 

{Q\P(x±)\Q,) cx (n\u(x±)\Cl) 
and the corresponding correlation function 

(n\T[p( Xl _)p(o)]\fi) cx (n\T[u(x±)u(o)}\n), 

which in turn yields the static quark-antiquark interaction energy . Up to an irrelevant factor we 
have after rotation to the Euclidean (r = 

cxp-t-LU (r)} = (fi|T [u(xf)u(0)] |Q), 

i.e. the static quark-antiquark potential is given directly by (the a = 6 = 3, [i = v = 2> component 
of) the vacuum polarization tensor 11"^, and not by the zero mass propagator with corresponding 
self-energy insertions as obtained in the standard Gaussian treatment. 



Similarly for evaluation of the adjoint Polyakov loop correlator (cf. Eq.(60)) we introduce the 
composite field v(x±) 

v(x ± )= f dx 3 $Ux)$>*{x). (73) 
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Figure 2: Ultralocality of the Polyakov loops. Wavy Lines denote Polyakov 
loop propagators, curly lines charged gluon propagators. 



Under center reflections v is invariant 



v(x±) -> v(x±), 



and expectation value and correlator of P & d(x±) are given by 



(n\p ad \n) ~ (n\ v \n) , (n\T[p ad {x ± )p ad (o)} \n) ~ («|t [«(a:±)»(o)] |n> 



(74) 



The system described by the effective action (69) exhibits remarkable properties already at the 



perturbative level. Most importantly the center symmetry is realized in the perturbative vacuum, 
i.e. in the ground state obtained by dropping all the terms containing the coupling constant g. Geo- 
metrical mass (Eq.(|7(i|)) and Aharonov-Bohm flux (Eq.(|7l"l)) are not affected by such a perturbative 
treatment. The perturbative ground-state is even under charge conjugation and the expectation value 
of the Polyakov loop vanishes 

{n pt \p( x± )\n pt ) = o, 

indicating an infinite free energy of a static quark. Indeed a perturbative analysis of the correlation 
function confirms this expectation. The novel aspects of perturbation theory in the center symmetric 
phase are related to the ultralocality property of the Polyakov loops. As a consequence of ultralocality, 
the Polyakov loop correlator is given by a one particle irreducible 2-point function, i.e. by the vacuum 
polarization rather than by the one particle reducible Green function of standard Gaussian variables. 
This is illustrated in Fig. |^. Thus the interaction energy of static charges is determined by the color 
neutral, jj, = v = 3 component of the vacuum polarization tensor 



V(r) = - -In 
L 



,rp n 33 (i Po>Pl> p 2 )e i K 



(27T) 



In one loop approximation (Fig. H) II is given by 



n 33 ( P ) = ^£<,3 

13 



d 3 q 



(2tt) 3 S 



D a Jl(q,q 3 )D ab ^(p 



g,qa) 



with the charged gluon propagator 



5 ab 



P 



pi 



M 2 + : 



P^Pv 



p\ 



M 



V2 



2tt 

a? 3, p 3 = — I n + - 



Using dimensional regularization in the evaluation of II 33 , ultraviolet divergencies occur only in 
the final sum over Matsubara wave- numbers. They are most conveniently regularized by a heat 
kernel method. Due to the antiperiodic boundary conditions of the charged gluon fields no infrared 
divergencies occur. The final result for II 33 yields the following asymptotic behavior of V 



V(r)«0ptr + T ln(r) + O(l) 



r/L > 1 
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with the perturbative value of the string constant 



2 / „ 7T 2 

^ = iV M2+ ^- 

The leading term in V is independent of the coupling constant. It is determined kinematically. At 
large separations, the dominant contributions to V arise from the singularities of the polarization 
tensor (cf. Fig. ||) closest to the real axis. They describe production of 2 charged gluons of mass 
M and minimal momentum (tt/L). The perturbative string tension is up to the factor l/L the 
threshold energy for this production process. With this argument we immediately can deduce the 
interaction energy of adjoint charges determined by the correlation function of the composite fields 
v(x±) (cf. Eqs.(||j74|)). Once more the large distance behavior of 14a is dominated by the threshold 



of two (charged) gluon production; in this case however also the ground state contributes. These 
considerations determine the asymptotic behavior of V^d up to a constant 



e -MUr) = (fi pt |P ad |ft pt ) 2 + 4 e ~ 2ffptr 

and thus 



Kd(r) ~ -e~ 2 ^ r . 
r 

The characteristic difference in the asymptotic behavior of the interaction energies of fundamental and 
adjoint charges respectively is due to the ultralocality and ultimately has its origin in the compactness 
of the Polyakov loop variables. 

As a final application, we discuss the effect of dynamical quarks on the interaction energy of static 
fundamental charges. Dynamical quarks give rise to an additional contribution to the order parameter 
field u(xj_) in Eq.@) 

f L r <-> 
u(x±) = J dx 3 (x) d 3 $ M (x) - ip{x) 7^73^ (z)} ■ 

The anti-periodic boundary conditions of the quark fields have been changed in the field redefinition 
Eq.(||) to 

^i,s 3 = L) = e- ivT3 / 2 i>(x ± ,x 3 = 0). (75) 

As for the charged gluons, this change in boundary conditions accounts for the interaction with the 
Aharonov-Bohm fluxes generated by the ultralocal Polyakov loop variables. As is clear from our above 
discussion of the interaction energy of adjoint charges, unlimited increase in the interaction energy 
with increasing separation of the static charges will avoided only if u (x±) develops a non- vanishing 
vacuum expectation value 

exp{-LV(r)}« \(n pt \u(0)\n pt }\ 2 . 
Applying again perturbation theory we have 

■jr(fipt|u (ll) | Opt) = -(ftpt|V> (x) y73^ 0*0 l°pt> = J J 7^2^" 

The divergent expression for the vacuum expectation value can be evaluated by performing the 
integral over transverse momenta in dimensional regularization and yields the result 



(Opt \u (x±) Opt) = — 7T- K 2 {nmL) 

n=l 



The parameter a characterizes the boundary conditions; a — 0, 1/2 denote periodic and anti-periodic 
boundary conditions respectively. For these values the order parameter u does not develop a vacuum 
expectation value. The value a = 1/4 denotes the quasi-periodic b.c. of Eq.(f75|) and yields a non- 
vanishing expectation value. For quark masses mL 3> 1 



(f2 pt |w (xjj |Q 



pt/ 



e 



-mL 



32 



and thus for asymptotic separations, the interaction energy of two static charges reaches the constant 
value 

V(r) « 2m . (76) 

It is remarkable that this correct form of the asymptotics of V appears within perturbation theory. 
In standard perturbation theory quark loops yield as in QED an Uehling type correction f3S|| to the 
Coulomb interaction SV cx g exp (— 2mr)/r. Perturbation theory in the center symmetric phase 
on the other hand yields the coupling constant independent result (|76|). Thus the mechanism of 
string breaking by pair formation is apparently present in the center symmetric phase already at 
the perturbative level. The calculation also displays the important role of the modification of the 
boundary conditions, i.e. the role of the Aharonov-Bohm fluxes. The string breaking mechanism 
would not be present if dynamical quarks satisfied standard anti-periodic boundary conditions, nor 
would it arise if the boundary conditions would not differentiate between the two color states of the 
quarks (t 3 in the b.c. ([75])). 



4.8 Polyakov Loops in the Plasma Phase 

In the last two sections we have identified the role of the Polyakov loop variables in the confined 
phase of QCD. We have displayed the relevance of these variables for both a phenomenological 
characterization of the confined phase and for development of a modified perturbative approach. 
We now will extend our investigations and will exhibit the particular role of these variables in the 
deconfined phase. The perturbative center symmetric phase with its signatures of confinement cannot 
be relevant for QCD at extensions smaller than L c = 1/T C or temperatures beyond the critical 
temperature T c . Not only do we expect the center symmetry to be broken at small extensions 
but also dimensional reduction to QCD2+1 to happen. Aharonov-Bohm fluxes induce anti-periodic 
boundary conditions and therefore yield a decoupling of the charged gluons if dimensional reduction 
takes place in the center symmetric phase (high temperature confining phase). The small extension 
or high temperature limit of the center symmetric phase is therefore QED2+1. In comparison to 
the dimensional reduction to QCD2+1, the resulting values for energy density and pressure would be 
reduced for SU(iV) Yang-Mills theories by a factor 1/(N + 1). This is in conflict with the results 
from lattice calculations. In order to reach the correct high temperature phase, the deconfinement 
phase-transition arising when compressing the QCD vacuum, must be accompanied by screening of the 
Aharonov-Bohm fluxes and simultaneously by a weakening of the increase in the geometrical mass. In 
the following we shall treat the Aharonov-Bohm fluxes (Polyakov loop variables) as phenomenological 
space-time independent quantities and as is appropriate for the deconfined phase, shift back the origin 
of these variables (cf. Eq.Q)) 

X = gLa 3 + ir. 

For simplicity we assume the geometrical mass M (cf. Eq.(|7(])) to vanish in the deconfined phase. 

The strength x is limited by the requirement of thermodynamic stability. By covariance, positive 
Casimir energy e at finite extension corresponds to a negative pressure of the corresponding system 
at finite temperature. The expression for the Casimir energy of the charged gluons 

71=1 

yields the following requirement for thermodynamic stability 

X < 1.51. (77) 

Complete screening (x = 0) of the Aharonov-Bohm fluxes, compatible with thermodynamic stability, 
is unlikely to take place in the deconfined phase. Such a system will exhibit, like perturbative QCD 



at zero temperature, a magnetic instability 19 . Spontaneous formation of magnetic fields on the 
other hand is prevented in the presence of sufficiently strong Aharonov-Bohm fluxes. Calculation of 
the Casimir energy for gluons subject to a color-magnetic homogeneous background field 

^ulbg = S a % lX2 H . 
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leads, for given extension L or temperature T to a lower limit for \. Figure [| displays the region 
of thermodynamic and magnetic stability. Obviously, due to the requirement of magnetic stability, 
the Stefan-Boltzmann limit (corresponding to x = 0) cannot be reached for any finite temperature. 
Identification of the Aharonov-Bohm flux with the minimal allowed values sets upper limits to energy 
density and pressure. They are shown in Figure ^ and are reminiscent of lattice data |39|| in their 
slow, logarithmic approach to the Stefan-Boltzmann limit. The finite value of the Aharonov-Bohm 
flux accounts for interactions present in the deconfincd phase and specifically gives rise to values of 
the interaction measure e — 3P which are in qualitative agreement with lattice data. Other quantities 
like the Debye screening mass are also affected by the non-vanishing Aharonov-Bohm flux. A one 
loop calculation 

m li = JT B * {£) > B 2 (x) = 1/6 - x(l - x) 

reproduces qualitatively the approximate linear dependence of m^i over a large temperature regime as 
observed in lattice calculations. Finally these results can also be used to estimate the change in energy 
density — Ae across the phase transition. In this phenomenological treatment the phase transition is 
accompanied by a change in strength of the Aharonov-Bohm flux from the center symmetric value it 
to a value in the stability region of Fig. H. The lower bound is determined by Eq. {f77 



-Ae = e(L c , X = tt) - e(L c , X = 1.51) < ^-jj . 

For deriving the upper bound the critical temperature must be specified. With T c « 270 MeV we 
obtain 

0.38 -L<-Ae< 0.53 



The lattice result 40 is within these limits 



Ae = -0.45-1 . 

It appears that to large extent, the dynamics in the deconfined phase is summarized by the dynamics 
of the Polyakov loop variables. The picture of increasing screening of the Aharonov-Bohm fluxes with 
increasing temperature seems to catch the essential physics of the thermodynamic quantities. It is 
remarkable that the requirement of magnetic stability which prohibits complete screening seems to 
determine the temperature dependence of the Aharonov-Bohm fluxes and thereby to simulate the 
non-perturbative dynamics in a semiquantitative way. 



4.9 Monopole Dynamics 

In the last sections we have investigated the dynamics of the Polyakov loop variables. We have de- 
scribed phenomenological implications and have demonstrated the important role of the compactness 
of these variables in the structure of the theory. In particular, the compactness offers the possibility 
of a perturbative approach to certain confinement related phenomena of QCD. For instance, when 
neglecting couplings to other gluonic degrees of freedom, the dynamics of the Polyakov loops is deter- 
mined completely by their compact nature and yields a gap in the spectrum which increases linearly 
with the extension of the system. This in turn implies a linearly rising interaction energy of static 
charges with a value of the string tension which becomes infinite in the continuum limit. If in a 
second step, gluonic couplings are taken into account perturbatively, the linear rise in the interac- 
tion energy persists with a value of the string tension which decreases with increasing extension of 
the system. Other qualitative features are properly described, such as the absence of confinement 
for adjoint charges or the string breaking mechanism induced by dynamical quarks. The formalism 
however also indicates that realistic values of the string tension are not attainable in this modified 
perturbative approach. A gap in the spectrum which rises linearly has to be present in sectors ex- 
cited by other gluonic degrees of freedom such as the order parameter field u(x±) defined in ([72]). 
Compactness of the Polyakov loop variables is not sufficient; variables like u(x±) also have to display 
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Figure 3: Regions of stability and instability 
in the (L, x) plane. To the right of the cir- 
cles, thermodynamic instability; above the 
solid line, magnetic instability. 




Figure 4: Energy density and pressure nor- 
malized to Stefan-Boltzmann values vs. tem- 
perature in units of Ams- 



similar properties. Origin or mechanisms for making such variables compact have not been identified. 
In the present axial gauge representation it is tempting to connect such mechanisms to the presence 
of monopoles whose existence is intimately linked to the compactness of the Polyakov loop variables. 
As we have seen, due to the quantization of the electric flux, the compactness leads to a magnetic 
field dominated ground state and exhibits the dual Meissner effect. All these characteristic properties 
of the QCD vacuum are however only realized as far as the Polyakov loop variables are concerned. 
Compactness of the Polyakov loop variables on the other hand implies the presence of singular field 
configurations whenever the Polyakov loop variables reach the border of their domain. Condensation 
of the points where the Polyakov loop passes through the center of the group is made likely by the 
compactness of these variables. The processes which control the appearance of these singular fields 
are not understood. Unlike 't Hooft-Polyakov monopoles of the Higgs model or instantons in QCD, 
the singular field configurations emerging in the gauge fixing process, in general, are not solutions to 
the classical field equations and the associated action can be arbitrary small. Their appearance is 
therefore not suppressed by a necessarily large contribution to the action. On the other hand, the 
entropy generated by the quantum fluctuations will be reduced when monopoles are present. Only if 
the fluctuations satisfy certain boundary conditions (cf. Eq.(|55|)) the action will remain finite. This 
interplay between monopoles and fluctuations presents a difficult theoretical problem. After decom- 
posing a certain configuration into singular and fluctuating components the action contains terms 
which are linear in the fluctuations. Furthermore a canonical description of this dynamics has not 
been derived. In our discussion of the monopole dynamics we will therefore adopt a phcnomenologi- 
cal point of view and on the basis of well understood properties of QCD we will attempt an indirect 
characterization of the dynamics. 

In the first part we will present arguments concerning the frequency of occurrence of singular field 
configurations in the partition function. To this end we make use of a relation between singular field 
configurations and instantons. As is easily verified, the Polyakov loop of a single instanton of size p 
(p <C L) is given by 
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Here and throughout this section we assume space-time to be Euclidean and use vector notation for 
x± and corresponding quantities. The above equation shows that, in axial gauge, a single instanton 
contains a north (P = 1) and south pole singularity (P — ~ 1) at its center and at infinity respectively. 
More generally it has be shown ^1, [l2| ^| that the topological charge v of a field configuration is 
given by the difference of the net northern and southern charge 



This relation implies, that field configurations which carry topological charge are necessarily singular. 
This is plausible. In the axial gauge one removes almost completely one space time component of 
the gauge field. The remaining 3 components (together with some "zero modes"), if continuous, are 
not sufficient to generate a non-vanishing topological charge. Thus in order to cover topologically 
non-trivial sectors, the ensemble of gauge fixed field configurations must contain singular fields. 
Furthermore, the action of singular fields with a topological charge, containing for instance a north 
pole and a south antipole, cannot be smaller than the instanton action ^r- The above relation can 
be used to estimate the density of monopoles um on the basis of the density nj obtained in either 
the instanton liquid model [Q or from lattice QCD jf6|. In order to relate instanton and monopole 
densities, we treat the instantons as independent; in this case the Polyakov loop is given by a product 
over the Polyakov loops associated with single instantons (or anti- instantons). For small instanton 
sizes /)<L, the number of instantons plus anti-instantons integrated over in the time integral of the 
Polyakov loop for fixed is of the order of Ln\p 3 . On the average, the number of instantons and 
anti-instantons will be the same and thus the phase accumulated in the time integral will be given by 
the fluctuations and is thus expected to be of the order of ny/ Ln\p 3 . When changing the position xj_ 
by an instanton size of p we expect, in the absence of correlations between the instantons, a different 
value of the phase which however is of the same order of magnitude. Thus, in this change of xj_ the 
Polyakov loop will typically pass \J Lnip 3 times through the center of the group which implies the 
following estimate for the monopole density (counting both poles and antipoles) at low temperatures 

n M oc (Lnip) 3 ^ 2 , (78) 
while for high temperatures, the integral involves at most one instanton and we therefore expect 

tim oc Ln\ , p > L. (79) 



Surprisingly, the result (78) implies an infinite monopole density at infinite extension or zero tem- 
perature. Although derived on the basis of a finite value of the instanton density, this result is to 
a large extent independent of the particular model; the characteristic L 3 / 2 dependence is mostly a 
consequence of the integration over x% which becomes ill-defined at infinite extension. It thus appears 
that irrespective of the temperature, axial gauge monopoles occur with non- vanishing or even infinite 
density. 

Beyond generation of monopoles via instantons, the system has the additional option of producing 
one type (north or south) of poles and corresponding antipoles only. No topological charge is asso- 
ciated with such singular fields and their occurrence is not limited by a bound on the action. Thus 
entropy favors production of such configurations their suppression could only come through their 
coupling to the quantum fluctuations. Since it appears that instantons may not be able to account 
for confinement |Q or, more precisely, for a realistic value of the string constant p6| , this additional 
option could be relevant. These entropy arguments also apply in the plasma phase. For purely 
kinematical reasons a decrease in the monopole density must be expected as the estimates ( |78| , [79| ) 
show. This decrease is counteracted by the enhanced probability to produce monopoles when, with 
decreasing L, the Polyakov loop approaches more and more the center of the group, as has has been 
discussed in the previous section (cf. Fig. [|). One certainly has to expect a finite density of singular 
fields also in the deconfined phase. In order for this to be compatible with the partially perturbative 
nature of the plasma phase and with dimensional reduction to QCD2+1 poles and antipoles may have 
to be strongly correlated with each other and to form effectively a gas of dipoles. 
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So far our qualitative arguments have dealt with the entropy associated with a finite monopole 
density. We now address dynamical issues. Singular fields may contribute only little to the action. 
For this reason we focus the discussion on the coupling of the quantum fluctuations 8 Ay, to the 



singular fields Sy (cf. Eq.(57)). We first consider the coupling of the charged quantum fluctuations 



to the Abelian monopole fields via the 4-gluon vertex, 

3 

8S= -g 2 j d 4 xY,^(x)-Stf>(x) a s (x ± ) • a s (x ± ). 

i—l 

At long wave-lengths the quantum fluctuations thus acquire a mass 8m by coupling to the Abelian 
magnetic field, 

<5m 2 =. 9 2 ij^a s (x ± ).a s (x ± ). 

To compute this mass term we assume the singular field a s (xj_) to be given by a superposition of the 
standard monopole fields of Eq. (B2 



N 

a s (x ± ) = Y a"' (x_l - x_u) 

i=l 



with positions of the monopoles and charges rii = ±1 and vanishing total magnetic charge, 

N 



i=l 



In this way the functional integral over the Polyakov loop variables in (|56|) has effectively been 
replaced by a summation over positions x_u and charges rii of these prescribed representative fields. 
Using standard identities from electrostatics, the mass term can be written as 



N N 



8m2 = Jv l, d " x T, l Y , "2", =~V £ ^•Ixx.-xx,-!. (80) 



4V Jy A —' Xj_ - Xj_i X_L - Xj_j V. 

i<j 

Similarly, neutral gauge field fluctuations acquire a mass Sm^ by coupling to the singular charged 
components, 

8mj=g 2 ± ^d 3 ^0 st (x ± ).0 s (x ± ). 

Evaluation of this expression by using the charged partner (Eq.Jsil)) of the Dirac monopole field 
yields 

8m\ = 8m 2 . 

Equation (^) displays the resistance of the system to proliferate production of monopoles for entropy 
reasons. The mass 8m 2 decreases with decreasing distance between monopoles of equal charge. 
Coupling to quantum fluctuations thus induces attraction between monopoles of opposite charges 
and repulsion between equal charge monopoles. If strong enough the attraction leads to monopole- 
antimonopole annihilation and ultimately to a vanishing 8m 2 . Computation of the mass requires 
knowledge of the correlations between monopoles. For an uncorrelated gas, the mass vanishes. If, 
for instance, the dynamics favors monopoles of opposite charge to be bound in dipoles and if the 
extension of the dipoles is much smaller than their separation, the expression (^0|) can be simplified 
and yields a gluon mass 

8m 2 = nyid 

determined by monopole density and average dipole size d. If applied to the deconfincd phase, it 
could provide a natural explanation for the appearance of a magnetic gluon mass. The temperature 
dependence of such a magnetic mass is determined by monopole density and dipole size. It is not 
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obvious that such a mechanism of mass generation is appropriate for the confined phase. Here, non- 
perturbative mechanisms must be operative which distinguish between the center-parity even and 
odd states (cf. Eqs.(59 6l])). Mass generation does not provide such a distinction, on the other hand, 
boundary conditions imposed by the presence of monopoles such as ( j55[) may do so. Fluctuations 
around singular fields have to satisfy these conditions otherwise infinite action results. Given a finite 
monopole density and assuming no particular correlations between the monopoles, this condition 
cannot be satisfied for long wave-length fluctuations, i.e. we expect fluctuations with wave number 

^ — ^min = n M 

to be dynamically suppressed. On the other hand, long wave length excitations associated with 
bilinears such as & (x)-<f>(x) which generate Z-even states are not necessarily affected by the boundary 
conditions and may be associated with excitation energies which remain finite for infinite extensions. 



5 Conclusions 

We have discussed QCD in the framework of gauge fixed formulations. We have provided a detailed 
analysis for both the Coulomb- and the axial gauge. The focus of our discussion has been on non- 
perturbative phenomena and we have emphasized the role of compact variables and of singular field 
configurations which both arise as a result of fixing the gauge. The appearance of such peculiar con- 
figurations is well understood and reflects obstructions in imposing these gauge conditions. However, 
the gauge dependence of these exceptional configurations remains disturbing and has been prohibitive 
for establishing the dynamical role of singular fields. Existence of obstructions in gauge fixing may 
be the only property which is common to the various gauge fixed formulations of QCD. It is easily 
seen that on a general level no statement concerning the importance i.e. the weight in the generating 
functional of such configurations is possible. One may for instance impose the Coulomb- or the axial 
gauge condition in the Georgi-Glashow model. These gauge conditions are independent of matter 
fields and so are all the geometrical properties of the fields associated with the particular gauge fixing 
procedure. Nevertheless, depending on the self-interaction of the Higgs particles, the system can be 
either in the Higgs or the confined phase. In the Higgs phase one does not expect these singular fields 
to be of any significance. The weight of the exceptional field configurations is thus determined by the 
dynamics. To clarify their relevance requires, dynamical studies beyond classification are required. 

In Coulomb- and axial gauge representation of QCD, the role of the exceptional fields is very dif- 
ferent. In Coulomb-gauge, Gribov horizons are given only implicitly and their presence affects all field 
configurations whose amplitudes exceed certain limits. In axial gauge only a particular set of degrees 
of freedom, the Polyakov loop variables, are subject to the compactness requirements. This makes it 
much easier to implement compactness in the axial gauge, where the Gribov horizon is determined by 
geometrical properties. On the other hand, compactness of variables other than the Polyakov loops 
has to emerge dynamically within the axial gauge. Our analysis has provided evidence that compact 
variables are the source of confinement- like phenomena. This is a result of approximative calculations 
carried out in both gauges. Neglecting, apart from the compactness requirement, all gluon self inter- 
actions the spectrum of the resulting Hamiltonians has been analyzed. In Coulomb-gauge, Gribov's 
approach amounts to replace the perturbative Hamiltonian of non-interacting massless gluons with 
an Hamiltonian of gluons with a momentum dependent effective mass which diverges in the infrared. 
In the axial gauge Hamiltonian, the compact Polyakov loop variables are treated like quantum me- 
chanical degrees of freedom moving in an infinite square well. Common to both formulations is the 
suppression of elementary excitations at long wave-length and a deviation from the cquipartition be- 
tween electric and magnetic field energy. They both give rise to interaction energies of static charges 
which increase with distance. In axial gauge, a linear rise is obtained with the string constant given 
by the lattice strong coupling limit. In Coulomb-gauge a complete analysis has not been carried out; 
plausibility arguments also point to a linear increase. The two approaches lead to distinctively dif- 
ferent descriptions of the ground states corresponding to the compact variables. The approximation 
by a momentum dependent gluon mass yields electric field dominance, the treatment of the compact 
variables as quantum mechanical degrees of freedom in an infinite square well gives rise to a magnetic 
vacuum. 
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Improvement in the description of the dynamics requires gluon self interactions and possibly 
coupling to dynamical quarks to be taken into account. Such a calculation has been described for 
the axial gauge, results in Coulomb-gauge do not exist. Confinement like properties get very much 
reduced in strength by coupling pcrturbatively the compact Polyakov loop variables to other gluonic 
variables. The linear rise in the interaction energy of static charges persists. However the string 
constant decreases with the extension of the compact direction. The QCD scale generated by the 
gluon loops does not at all affect the string constant. When dynamical quarks are coupled to the 
gauge fields, string-breaking occurs. It is surprising that the minimal non-perturbative dynamics 
incorporated into the compactness of the Polyakov loops is sufficient to lead to this non-perturbative 
mechanism. Softening or loss of confinement must be expected to occur also in Coulomb-gauge. The 
interaction energy receives loop corrections from transverse gluons corresponding to the screening 
contribution to the running coupling constant. These corrections are likely to lead to an interaction 
energy which remains finite for large distances since in Gribov's treatment the gluon propagator 
actually vanishes in the infrared. Perturbative treatment of gluonic self-interactions appears to be 
inappropriate for improving the description of the confinement like phenomena linked to the presence 
of compact variables. 

Gauge fixed formulations are not only important for developing approximation schemes beyond 
standard perturbation theory, they also provide a useful framework for more phenomenologically 
oriented analyses of QCD. In particular axial gauge offers the possibility for a phenomenological 
approach to the Polyakov loop dynamics. These variables serve as order parameter fields to charac- 
terize the phases of QCD via the realization of the center symmetry. In the confined phase with the 
center symmetry realized, the states of the system can be classified (for SU(2)) according to their 
behavior under center reflections. Hadronic states like glueballs have even center parity. Excitation 
energies of states of negative center parity, with vanishing momentum in the compact direction, in- 
crease with the extension of the system. These two sectors of Yang-Mills theories are completely 
decoupled below the phase transition and merge into an interacting system above. Thus the phase 
transition cannot be understood as a reordering in the glueball sector, it involves a sector of the 
Hilbert space which is invisible above the critical extension or below the critical temperature. The 
Polyakov loops also characterize to some extent the high temperature phase. Analysis of thermody- 
namic and magnet stability limits the range of values the Polyakov loop can assume in the phase of 
broken center symmetry which in turn imply limits on the deviations of energy density and pressure 
from the Stcfan-Boltzmann values. 

The most severe limitations in further developing gauge fixed formulation are related to the pres- 
ence of singular field configurations. Unlike instantons or the 't Hooft-Polyakov monopole, singular 
field configurations in general do not satisfy the classical field equations. The action associated with 
such fields may be arbitrarily small. Semiclassical treatments are not directly applicable. Singular 
fields and fluctuations can mix. In general fluctuations around a singular field destroy the delicate 
balance of Abelian and non-Abelian contributions to the field strength resulting in an infinite action. 
There are no standard tools available to handle the difficult constraints imposed on the fluctuations 
by finiteness of the action. Essentially nothing is known about the treatment of singular fields in 
the canonical formalism. It is not obvious how to include them in an Ansatz for the QCD ground 
state wave-functional. Simple QED like structures for the wave-functional would exclude such fields 
by erroneously assigning to them an infinite energy. It seems unlikely that the issue of singular fields 
can be avoided. Smooth, but topologically non-trivial configurations must be in general expected to 
give rise to singular fields after gauge fixing. Once more one may approach these difficult dynamical 
issues also from a phenomenological point of view. We have demonstrated this procedure by inferring 
from the high temperature limit of Yang-Mills theories, from results of lattice QCD or of models like 
the instanton liquid model some basic properties of singular fields such as the their density, and its 
temperature dependence or their dynamical correlations. Despite some insights gained in this way 
we are still far away from a sufficiently complete phenomenological characterization of the dynamics 
of singular fields. 
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